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Abstract 

We derive the parity- and time-reversal- violating nuclear interactions stemming from the QCD 9 
term and quark/gluon operators of effective dimension 6: quark electric dipole moments, quark and 
gluon chromo-electric dipole moments, and two four-quark operators. We work in the framework 
of two-flavor chiral perturbation theory, where a systematic expansion is possible. The different 
chiral-transformation properties of the sources of time-reversal violation lead to different hadronic 
interactions. For all sources considered the leading-order potential involves known one-pion exchange, 
but its specific form and the relative importance of short-range interactions depend on the source. 
For the 6 term, the leading potential is solely given by one-pion exchange, which does not contribute 
to the deuteron electric dipole moment. In subleading order, a new two-pion-exchange potential 
is obtained. Its short-range component is indistinguishable from one of two undetermined contact 
interactions that appear at the same order and represent effects of heavier mesons and other short- 
range QCD dynamics. One-pion-exchange corrections at this order are discussed as well. 



1 Introduction 



There are three sources of CP violation in the lowest-dimension operators of the Standard Model. Due 
to its flavor-changing properties, the phase of the CKM matrix is the best known and investigated 
source (see, for example, Ref. [T]). The recent discovery of neutrino masses opens up the possibility of 
analogous leptonic CP violation [2] . The third possible source comes from the 9 term in QCD [3] . Flavor- 
conserving quantities, in particular hadronic, nuclear, and atomic electric dipole moments (EDMs), are 
sensitive probes of this qualitatively different mechanism of CP violation. EDMs require the simultaneous 
violation of parity (P) and time-reversal (T); for a review, see Ref. [3]. 

Years of careful experimental investigation have set stringent bounds on the neutron EDM, \d n \ < 
2.9 • 10 -26 e cm [5]. This upper limit points to a tiny value for the 9 parameter, (9<10~ 10 , despite 
the naive-dimensional-analysis expectation that it be of 0(1). Besides the neutron EDM, there are 
measurements on atomic EDMs. The leading contribution in a paramagnetic system comes from the EDM 
of the unpaired electron. In contrast, a diamagnetic atom has zero total electron angular momentum, 
and the main contribution to the atomic EDM comes from the nuclear Schiff moment (SM) , the residual 
electron-nucleus interaction generated by the distribution of the EDM throughout the nuclear region. A 
less stringent bound on the proton EDM, \d p \ < 7.9 • 10 -25 e cm, is extracted from the EDM of the 199 Hg 
atom [5] through a calculation of the nuclear SM [7J. Atomic EDMs can also receive contributions from 
higher nuclear moments that violate T, such as the magnetic quadrupole moment (MQM). 

The unnaturally small value of 9 leaves room for other flavor-diagonal sources of T violation, which 
have their origin in physics beyond the Standard Model, at a high-energy scale Mt- Well below the new 
physics scale, these effects manifest themselves in interactions between Standard Model fields represented 
by higher-dimension T- violating (TV) operators, suppressed by more and more powers of Mp. We expect 
the most important TV effects to be captured by the operators of lowest dimension [5] |H1 QUI H] , the 
dimension-6 quark electric and chromo-electric dipole moments (qEDM and qCEDM), gluon chromo- 
electric dipole moment (gCEDM), and TV four-quark (FQ) operators. 

A recent renewal of the longstanding interest in EDMs comes from a new generation of experiments, 
which will improve the precision of EDM observables significantly. It is expected that the current bound 
on the neutron EDM will be pushed down to 10~ 27 — 10 -28 e cm [I2 j using high-density ultra-cold neutron 
sources at SNS Q2] and ILL+PSI [2]. It has also been proposed that the EDM of charged particles could 
be measured in storage- ring experiments [15]. In addition to the proton, we could see a measurement 
of the deuteron EDM with a projected precision of \dd\ ^ 10 -29 e cm. A measurement on 3 He is also a 
possibility. These measurements probe dimension-6 sources at scales comparable to the scale reached at 
the LHC, so they could very well turn up positive signals. 

Calculating hadronic and nuclear EDMs (and higher moments) directly from QCD with TV sources is 
a daunting task, although progress has been made for nucleons, in the case of 9, using lattice simulations 
[16] . An alternative is to use a low-energy effective field theory (EFT) of QCD, chiral perturbation theory 
(ChPT) [TjJ [18] (for a review, see for example Ref. 19 ]). In this case, we can describe hadronic and 
nuclear observables in a controlled expansion with minimal dynamical assumptions, where the symmetries 
of QCD — in particular the chiral SUl(2) x SUr(2) symmetry for two quark flavors — are respected 
order by order. While the simplest T-violating (TV) interactions have been known for a long time 
[20] [21] [22] [23] [24] , a comprehensive analysis — similar to the well-known T-conserving (TC) sector — of 
the TV chiral Lagrangian from the 9 term [25j and the dimension-6 sources [26j has only been performed 
recently. 

In the case of the nucleon EDM, one finds that, for T violation from the 9 term, the leading contri- 
butions come from the pion cloud, where the pion couples to the nucleon via a non-derivative P- and 
T-odd interaction, and from shorter-range interactions. The former is purely isovector, and provides an 
estimate of the EDM: the characteristic lnm^, with the pion mass, is not expected to be canceled 
by short-range contributions and thus provides the bound on 9 mentioned above [20] [21] [23] . One can 
extend the calculation to the full electric dipole form factor (EDFF) as well [23], where in leading orders 
the radius, a contribution to the nucleon SM, is predicted [35] ■ However, because of two possible short- 
range structures, isoscalar and isovector, even measurements of both neutron and proton EDMs would 
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be insufficient to determine the three unknown parameters that appear at ieading order. 

When dimension-6 sources are included in the picture, one finds that the qCEDM generates a nucleon 
EDFF that is not distinguishable from the one stemming from 6 24 . Furthermore, while the momentum 
dependence of the EDFF is qualitatively different if the qEDM or the gCEDM, rather than 9 or the 
qCEDM, are the dominant sources of T violation, a measurement of the neutron and proton EDMs 
alone can be fitted equally well by any of these sources. It is clear, then, that to pinpoint the dominant 
mechanism(s) of T violation at high energy, more observables are needed, and nuclear EDMs are natural 
candidates. 

Nuclear EDMs and other moments receive various contributions. There are, of course, contributions 
from the individual nucleons' EDMs. In the deuteron the isovector component cancels, while in 3 He 
one can expect a cancellation between the contributions of the two nearly anti-aligned protons. Thus, 
nuclear EDMs, in particular the deuteron's [27J 12H]) are sensitive to a different combination of hadronic 
TV parameters than the neutron EDM. However, these one-nucleon contributions are modified from their 
"in-vacuum" counterparts because the nucleons are not free but bound in the nucleus. There are many- 
nucleon effects that are TV. First, the TV component of the pion cloud can generate a TV pion-exchange 
interaction among nucleons, and no symmetry forbids interactions of shorter range, either. These TV 
nuclear forces will mix-in components of the nuclear wave function that do not appear in the absence 
of TV. It is a source of polarization effects for the entire nucleus. Second, there may be multi-nucleon 
contributions to the TV coupling of the photon; such TV currents can be generated by either pion 
exchange or shorter-range dynamics. 

TV one-pion exchange (OPE) has long been recognized as an important component of the TV two- 
nucleon (NN) potential, and expressed [2^1 [30] m terms of three non-derivative pion-nucleon couplings 
[3T] , associated with isospin 1 = 0, 1,2. So far, the analysis of TV nuclear effects has been based on 
tree-level potentials where OPE is sometimes supplemented by the single exchange of heavier mesons, 
the eta [32], rho [33], and omega [33] being most popular. Allowing sufficiently many couplings of these 
mesons to nucleons one can produce [34] the most general short-range TV NN local interaction with one 
derivative [35 . This is the TV analog of the DDH approach [36] for nuclear TC P violation (PV). 

The contributions from such potentials to the deuteron and 3 He EDMs have been calculated in the 
literature under various assumptions. OPE from the 1 = 2 TV pion-nucleon coupling does not contribute 
to the NN system at tree level. It was noticed early on [37] that OPE from the 1 = TV pion-nucleon 
coupling does not contribute to the deuteron EDM, either, but it does for 3 He, where it was estimated 
with a phenomenological strong- interaction potential [38] . The deuteron EDM that arises from an / = 1 
TV exchange of either pion- or shorter-range, together with a separable strong-interaction potential, was 
calculated in Ref. [37]. The effects of OPE on the deuteron EDM and MQM were calculated using both 
zero-range and phenomenological strong-interaction potentials in Ref. [39]. More recent calculations of 
the deuteron EDM and MQM [34] and of the 3 He EDM [40] have considered other TV contributions 
besides those from the TV potential, and used more modern, "realistic" strong-interaction potentials. 
Meson-exchange currents were found small in the deuteron [34], and neglected in 3 He [40]. The TV- 
potential contributions are consistent with earlier results; they are dominated by OPE from the 1 = 1 
pion-nucleon coupling in the case of the deuteron [34], and from all three pion-nucleon couplings in the 
case of 3 He [40] . 

TV moments of heavier nuclei are more difficult to calculate. It has been argued [29 that the TV 
potential can lead to an enhancement over the nucleon EDM thanks to the near-degeneracy of levels of 
opposite parity, while meson-exchange currents are comparatively small. The size of the effect can be 
estimated through the single-particle potential obtained by averaging the NN potential over a closed 
nuclear core. The OPE from the I = 0,2 pion-nucleon couplings are proportional to the nuclear J3, 
(N — Z)/A [29] [301 141], while the OPE from the 1 = 1 coupling does not have such a suppressing factor 
[30] |4T] . It has also been found that the matrix elements from the rho and omega are small compared to 
the I = 1 pion contribution 33 . EDMs and MQMs (for example from 1 = OPE [39]) and SMs g3] of 
several interesting nuclei have been estimated. A sample of recent SM calculations can be found in Refs. 

Bug. 

There are, of course, other nuclear tests of TV, see for example Ref. [44]. The most promising for 
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effects of the TV nuclear interaction seems to be neutron scattering [35] . On the proton [33] and deuteron 
[4"T] . TV neutron scattering is again dominated by OPE, but sensitive mostly to the / = 0, 2 and I = 0, 1 
couplings, respectively. For heavy nuclei, one can again obtain estimates using the single-particle potential 
[301144]. 

For consistency, we would like to describe nuclear TV observables in the same framework used for the 
calculation of the nucleon EDM. The non-analytic behavior of the nucleon EDM in and the domi- 
nance of OPE in nuclear observables point to the need of a framework that can account for both effects 
simultaneously, with chiral symmetry playing a central role. In fact, different sources of T violation have 
different transformation properties under chiral symmetry. As a consequence, the relative importance of 
various pion-nucleon and short-range interactions is not the same for all sources. Here we use the chiral 
Lagrangian built in Refs. [25, 26], where TV interactions stemming from the 6 term and the dimension-6 
sources were constructed and ordered according to the same power counting used to order TC interactions 
in ChPT. 

In addition to consistency between one- and few- nucleon TV interactions, nuclear TV also requires 
consistency between TV and TC forces, in order there to be no mismatch in the off-shell behaviors of 
the various ingredients. Of course, off-shell effects are dependent on the choice of fields, while physical 
quantities are not, provided the same choice of fields has been made throughout the calculation. As far as 
TV nuclear interactions are concerned, phenomenological TC models bring additional uncertainties, such 
as the choice of zero-range or finite-range interactions and the role of heavy mesons. On the other hand, 
ChPT has been extended to multi-nucleon systems [IS], leading to the derivation of TC nuclear forces 
and currents. This opens up the possibility of describing all necessary ingredients in a single framework. 

The goal of this paper is to provide the first step in the extension of TV interactions in the EFT to the 
multi-nucleon sector where OPE is treated non-perturbatively. Some of us have recently looked at the 
deuteron's PV, TV electromagnetic moments in an EFT where pions are treated in perturbation theory 
[28] , In this case, the size of uncertainties is set by the relatively low scale where OPE becomes significant. 
Treating OPE non-perturbatively extends the EFT to higher momenta and improves convergence. The 
TV nuclear potential is the most important ingredient in this extension, and here we derive it for the 
most important TV sources. This is the TV potential to be used, for example, with the TC, parity- 
conserving (PC) potentials from Refs. [49] [50] HU [52] [53] [54] [55]. The construction here is similar to 
that of the TC PV potential [H [57] [58], which extends the EFT from the TC, PV one-nucleon sector 
[59] to multi-nucleon systems. Such a framework provides an alternative to the DDH approach [36], 
allowing for a mo del- independent analysis of nuclear TC P V phenomena [60] [61] . Our present TV PV 
EFT framework stands in respect to previous approaches like this TC PV EFT framework with respect 
to the DDH approach. 

As we are going to see, the form of the TV potential at a given order in the chiral expansion depends 
on the source. For the dimension-6 sources, all non- derivative pion-nucleon couplings appear at leading 
order, although these sources differ in the relative strength of the 1 — 2 pion-nucleon coupling and 
short-range interactions. For these sources, the leading-order potential is sufficient for most applications 
we envision. The situation is different for the 9 term: only the 1 = pion-nucleon coupling, which 
is suppressed in some cases of interest, appears at leading order. For this source, we thus derive the 
potential up to subleading order. We show that new elements appear with respect to phenomenological 
treatments, such as two-pion exchange (TPE) at the same level as short-range interactions representing 
heavier- meson exchange. This richness is a blessing, as it potentially reveals the TV source [24 ] 128 ] , 

The rest of the paper is organized as follows. In the next section, Sect. [5] we present a summary 
of basic ChPT ideas and give the TC and TV chiral Lagrangians needed in the following sections. For 
processes involving momenta below M nuc ~ 100 MeV, pion degrees of freedom can be integrated out 
and the dominant TV contact interactions (the leading-order TV potential in the so-called pionless EFT) 
are obtained in Sect. [3] The nuclear potential in ChPT, which should apply beyond M nuc , is then 
presented in momentum (Sect. [3]) and coordinate (Sect. [SJ spaces. (We relegate details of the Fourier 
transformation to the Appendix.) In Sect. [B]we discuss the size of different components of the potential 
and compare them with phenomenological forms. We draw our conclusions in Sect. [7] 
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2 Chiral Perturbation Theory 



QCD is characterized by an intrinsic mass scale Mqcd ~ 1 GeV. At momenta Q comparable to the 
pion mass, Q ~ m v <C Mqcd, interactions among nucleons and pions are described by the most gen- 
eral Lagrangian that involves these degrees of freedom and that has the same symmetries as QCD. A 
particularly important role at low energy is played by the approximate symmetry of QCD under the 
chiral group SUl(2) x SUr(2) ~ 50(4). Since it is not manifest in the spectrum, which instead exhibits 
an approximate isospin symmetry chiral symmetry must be spontaneously broken down to the isospin 
subgroup SUl+r(2) ~ 50(3) . The corresponding Goldstone bosons can be identified with the pions, 
which provide a non-linear realization of chiral symmetry. 

Chiral symmetry and its spontaneous breaking strongly constrain the form of the interactions among 
nucleons and pions. In particular, in the limit of vanishing quark masses and charges, when chiral 
symmetry is exact, pion interactions proceed through a covariant derivative, which in stereographic 
coordinates 7? for the pions is [52] 

£>„t? = D~ l d^ (1) 

with 

D = l + if 2 /F^ (2) 

and ~ 186 MeV the pion decay constant. One can also construct the covariant derivative of this 
covariant derivative, 

2 

VvD^TT = dyD^TT + — (TT A,7T • D^TT - D V W 7? • Dfjf) , (3) 
"it 

and so on. Nucleons are described by an isospin-1/2 field N, and we can define a nucleon covariant 
derivative 

V»N = (fy + JLf • (tt x D m tt)J N, (4) 

where Tj, i — 1, 2, 3, are the Pauli matrices in isospin space. We also define through NV^ = T>N, and 
use the shorthand notation 

v%_ ee v ± v n \ v>ijr ± = v^v" + v^v^ ± v^v v ± v^v (5) 

and 

ee t{D^ ± 2?t »Ti, t.{D^ ± V v ± ee t{D^V v + V^V^n ± V^t{D v ± V^t{D^. (6) 

When acting on a nucleon bilinear of non-zero isospin, for example T>^{NfN) 1 the covariant derivative 
is meant to be in the adjoint representation, that is, the isospin matrix in Eq. Q should be replaced by 

(V) ik = ie i i k . 

At Q ~ 771^ -C m^r, the nucleon mass, nucleons are essentially non-relativistic; as such the only 
coordinate with which their fields vary rapidly is v ■ x, where v is the nucleon velocity, — (1,0) in the 
nucleon rest frame. It is convenient therefore to use a heavy-nucleon field from which this fast variation 
has been removed |63j . This simplifies the gamma-matrix algebra, leaving only the spin operator 5^, 
where — (0, <?/2) in the nucleon rest frame. Below we use the subscript _L to denote the component 
of a four-vector perpendicular to the velocity, for example 

V lt± = V ll -v lt vV. (7) 

Baryon states above the nucleon, such as the delta isobar, can be included in EFT along similar lines, 
but for simplicity we do not include them here. 

The ChPT Lagrangian in the chiral limit includes all the interactions made out of D^if, N and their 
covariant derivatives that are chiral invariant. Chiral symmetry is explicitly broken, however, which 
introduces pion interactions that might not include derivatives, but are proportional to powers of the 
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symmetry-breaking parameters. Since the explicit breaking of chiral symmetry is small, chiral-symmetry- 
breaking operators can be systematically included as a perturbation on the chiral-invariant Lagrangian. 
Their forms are also not arbitrary, being instead determined by the chiral transformation properties of 
their progenitors in the QCD Lagrangian. The construction of chiral-symmetry-breaking operators from 
quark masses and "hard" electromagnetic interactions (those from photons with momenta beyond the 
EFT regime) is extensively treated in Refs. [62] EI]- "Soft" interactions via an explicit photon field 
appear in gauge-covariant derivatives, 



(D^TTa) -> — (d^Sab - eA M £ 3a6 ) 7Tfc, (8) 



1 

D 

-\ ■ * 1 + 7 " 3 

F 2 



®H + T^' x D ^) - ieA ^- 



N, (9) 



where e is the proton electric charge, and in gauge-invariant interactions built from the photon field 
strength 

= d^Ay - d v A^ (10) 

Approximate chiral symmetry, together with the heavy-baryon formalism, allows us to systematically 
expand observables in the mesonic and one-nucleon sectors in powers of Q/Mqcd, where Q is the typical 
momentum of the process under consideration. The ChPT Lagrangian contains an infinite number of 
terms, which can be organized using an integer "chiral index" A and the number / of fermion fields 

[171152]: 

oo 

£=££4 A) > (id 

A=0 / 

where A = d+ //2 — 2 > 0, with d the number of derivatives, powers of the pion mass or of the electric 
charge. For processes with at most one nucleon, A = 0, 1, all momenta and energies are typically ~ Q. 
The contribution of a diagram to the amplitude T can then be estimated by 

T oc Q v J 7 (Q/fx), (12) 

where T is a calculable function, /j, is the renormalization scale, and the counting index v is 

v = A-2C- A + 2L + J2 A >- (I 3 ) 

i 

Here, C — 1 and L are respectively the number of connected pieces and loops in the diagrams, and i 
counts the number of insertions of vertices from Ck ■ . From Eq. (|13[) it is apparent that diagrams with 
increasingly higher number of loops and non-vanishing-index interactions are increasingly suppressed, 
leading to a perturbative expansion. Assigning to loops a characteristic factor Q 2 /(Att) 2 and using naive 
dimensional analysis [T7 1 l65 | l9"] to estimate the EFT parameters, the suppression scale is Mqcd ~ 2-kF^. 
Note that in this sector of the theory nucleon recoil is a subleading effect: the nucleon is nearly static. 

The ChPT power counting formula (IT5|) cannot directly be applied to processes with A > 2 [55J [35] . 
Indeed, in diagrams in which the intermediate state consists purely of propagating nucleons — which are 
called "reducible" — the contour of integration for integrals over the 0th components of loop momenta 
cannot be deformed in way to avoid the poles of the nucleon propagators, thus picking up energies 
~ Q 2 /m,N from nucleon recoil, no longer a subleading effect, rather than ~ Q. There is also an extra 
factor of 4-7T. These diagrams are therefore enhanced by factors of Attitin/Q with respect to the ChPT 
power counting that assigns Q 2 /(An) 2 to a loop, and the need to resum them leads to the appearance of 
shallow bound states in systems with two or more nucleons, nuclei. Diagrams whose intermediate states 
contain interacting nucleons and pions — "irreducible" — do not suffer from this infrared enhancement, 
and in them nucleon recoil remains a small effect. Reducible diagrams are thus obtained by patching 
together irreducible diagrams with intermediate states consisting of A free-nucleon propagators. Calling 
V the sum of all irreducible diagrams, the amplitude can be written schematically as 

T = V + VGqV + VGqVGqV +... = V + VGqT, (14) 
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where Go is the free-nucleon, non-relativistic Green's function. Equation (|14[) is just the Lippmann- 
Schwinger equation, which is formally equivalent to a Schrodinger equation with a potential V. 

Naive dimensional analysis suggests [66] that irreducible diagrams follow the ChPT power counting 
rule (fT3|) with C > 1. While this is true for pion-exchange diagrams, the situation is more complicated 
for contact interactions. In fact, it can be shown that the iteration of the singular one-pion exchange 
requires for renormalization at the same order a finite number of / — 4 interactions, some of which are 
less suppressed than expected on the basis of naive dimensional analysis [671 168] . On the other hand, 
corrections, which should be perturbative, are expected to still conform to dimensional analysis [69l [70] , 
Since the TV potential is very small, it should be amenable to an expansion in powers of Q/Mqcdi 
with different contributions organized according to their chiral index v 1 or, equivalently, according to the 
number of inverse powers of Mqcd ■ 

In Sect. H]we compute the TV nuclear potential, in the case of the QCD 8 term up to order v = 3, 
which means up to 0(Q 2 /Mq CD ) with respect to the leading piece. Such a calculation requires the 
knowledge of the TC and TV ChPT Lagrangians up to A = 2 and A = 3, respectively. In the remainder 
of this section we present the relevant interactions. Throughout, we use nucleon field redefinitions to 
eliminate nucleon time derivatives from subleading interactions. 



2.1 T-Conserving Chiral Lagrangian 

The calculation of the TV potential in ChPT requires certain TC interactions with / = 0, 2, which we 
list here. (A more complete list can be found in the literature, for example Refs. [HJ EE1 HI] ) These 
interactions stem from the quark (color-gauged) kinetic and mass terms in the QCD Lagrangian. 
The leading TC chiral Lagrangian has chiral index A = and is given by 

4<2,t = ■ D ^ - ^ 2 + Niv ■ VN - ■ NtS^N, (15) 

Li Li j_J _t TT 

where gA is the pion-nucleon axial coupling, gA — 1.27. The pion mass term originates in explicit chiral- 
symmetry breaking by the average quark mass m = (m u + m c [)/2 and, by naive dimensional analysis, 
OimMQCD)- Neglecting for the moment isospin-breaking operators, at chiral order A = 1 the 



m 
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relevant Lagrangian consists of 

4< 2 ti = -^—NV 2 X N + (iv ■ Dtt) ■ NfS ■ V-N + Am N (l - NN. (16) 

} ^ A ' 1J 2m N x F^rriN \ F%D J 

Here the first two terms are the nucleon kinetic energy and a relativistic correction to the pion-nucleon 
coupling, the coefficients of both operators being fixed by Galilean invariance. The third term is the 
nucleon sigma term with a coefficient Amj/ = 0(to^/A/qc_d). At the next chiral order, A = 2, 



2g A 



2tt 2 \ 



■ NtS^N. (17) 



The first term is a correction to the pion mass, Am 2 — O (rri^/MQ CD ). The second term represents 
further relativistic corrections to the gA term in Eq. (|15l) . The constraints imposed by Lorentz invariance 
on Eqs. (fTo]) and (fT?)) agree with the results of Ref. [7T|, once a field redefinition is used to eliminate 
time derivatives acting on the nucleon field from the subleading A = 1 and A = 2 Lagrangians. The 
operator with coefficient ca = O (1/Mqcd) m Eq. (TT71) is a contribution to the square radius of the 
pion-nucleon form factor, while d,A = O (m 2 r /MQ CD j is a chiral-symmetry-breaking correction to gA [51j . 
which provide the so-called Goldberger-Treiman discrepancy. 

Isospin-breaking operators in the chiral Lagrangian |64] stem from the quark mass difference m c i — m u = 
2ms and from quark coupling to photons through the fine-structure constant a om = e 2 /47r. Here, 
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for simplicity, we count e ~ 1/3 as 0(1) and a em /47r as 0(m^/MQ CD ), since numerically a cm /47r ~ 
em 3 /(27rF' 7r ) 3 . Isospin- violating terms first contribute to the A = 1 Lagrangian, 

r W $ m l (-2 2\ , Sm N ,-y ( 2tt 3 _ A 



while at order A = 2, 



.(2) _ 5m 2 5m N - 
L f<2,Tf- 2D* 3+ 2 



T 3 + ^( 7r 3^-T-7f 2 T 3 ) 



V + ^ I D M 7T3 - ^7? • D P 7? ) V.S'" .V. 



2tt 3 



(19) 

Here 5m\ — (a cm MQ CD / Air} is the leading electromagnetic contribution to the pion mass split- 
ting, while the quark-mass-difference contribution, 5m\ — (e 2 rti^/MQ CD ) , is smaller by a power of 
em-ir /Mqcd- The pion mass splitting, m 2 ± — m 2 = <5r7i 2 — (5m 2 = (35.5 MeV) 2 [72], is dominated by the 
electromagnetic contribution. The nucleon mass splitting, m n — m p = 5m N + Sm N = 1.29 MeV [72j also 
receives contributions from electromagnetism and from the quark masses. In this case, the quark-mass 
contribution 5m n is expected to be the largest. By dimensional analysis 5m jy = 0{em\jMQCD\ and 
lattice simulations estimate it to be 5m^ = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV 73 , which is in agreement 
with an extraction from charge-symmetry breaking in the pn — > dir° reaction |74j . The electromagnetic 
contribution is 5m^ = {a cra MQCD/^), that is, (em^/MQ CD ) and about the 20% of 5mN- Using 
the Cottingham sum rule, 5m n = —(0.76 ± 0.30) MeV [75], which is consistent with dimensional analy- 
sis. The operator with coefficient Pi — O (eto^/Mq^J is an isospin- violating pion-nucleon coupling. At 
present there are only bounds on f3\ from isospin violation in nucleon-nucleon scattering. For example, a 
phase-shift analysis of two-nucleon data gives f3i = (0±9) • 10 -3 [51] [52], which is comparable to estimates 
of /3i from tt-tj mixing. 

For the solution of the Lippmann-Schwinger equation, it is convenient to eliminate the nucleon mass 
difference m n — m p from the nucleon propagator and from asymptotic states. This result can be accom- 
plished through a field redefinition, defined in Ref. [53]. After the field redefinition, Eqs. (TT8"]) and (TT9")) 
become 



f<2,Tj - ~2~D2 \ 5m * - §m N) (5 -^3) - ^TTjf 7 ^ - ( Sm N + 5m N ) (7? x v ■ Dtt) 3 

9ASmN ie 3ab n a Nr b S -V_N + ^- (D,n 3 - -^v? • D,n) NS^N. (20) 



£ (1,2) 



F,m N oa ° a " ' F v V F*D 

We will incorporate isospin-breaking effects in the potential using the Lagrangian (|20|) . 



2.2 T- Violating Chiral Lagrangian Prom 9 

The lowest-dimension TV operator that can be added to the TC QCD Lagrangian is the dimension-4 
9 term. With an appropriate choice of the quark fields q = (u,d) T , the 9 term can be expressed as a 
complex mass term |76j . 

£;T4 = m*9 qij 5 q, (21) 

where = m u md/ (m u + md) = (ui^/Mqcd) and 9 is the QCD vacuum angle, here already assumed 
to be small as indicated by the bound on the neutron EDM, 9 < 10~ 10 . The 9 term transforms under 
chiral symmetry as the fourth component of an 50(4) vector P — (qrq, qi^/sq), whose third component is 
responsible for quark-mass isospin violation [64] . TV from the 9 term and isospin violation from the quark 
mass difference are therefore intrinsically linked; this link appears in certain relations [20] [25] between 
the coefficients of TV and isospin-breaking operators in ChPT through a coefficient p = (1 — e 2 )9/2e. 

The pion-nucleon TV Lagrangian from the QCD 9 term was constructed in Ref. [25J. The leading 
/ < 2 TV interaction generated by the 9 term appears at A = 1 and consists of an isoscalar pion-nucleon 
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coupling, 

^fl 2 , n = -^-^N. (22) 

The relation to isospin-breaking operators implies that go can be expressed in terms of the quark-mass 
contribution to the nucleon mass difference, go = pSmpf = O (0tti^/Mqcd)- Increasing the chiral index 
by one, we find a single term, 

4=2,74 = -^*- D ^ NS^N, (23) 

where h is an isoscalar two-pion-nucleon coupling, which is related to the coefficient fi\ in Eq. (|20[) by 
ho = pPi = O ((9to^/Mq C £,). The A = 3 operators relevant to the calculation of the TV potential at 
order v = 3 are 

MS) 1 f 9i - Amjy 5ml \ (, 27?2 \ ^tat 

_ J_ / Ago 5K\ / _ . X / _ _2^\ . ^ 

+ 8^75 {* ' ^ + 2 ( x - 1) ^ ' ^ ^ ^ V »- N } • ^ 

where once again we eliminated operators containing nucleon time derivatives and neglected multi- 
pion operators. Equation (|24[) is obtained after rotating away the subleading pion tadpole operator 
— (po^w 2 r /25mjv)'7r3 from the mesonic TV Lagrangian, as detailed in Ref. |25j . The first operator in Eq. 
(|24|) is the most important contribution of the 9 term to the isospin-breaking TV non-derivative coupling 
n^NN. Its coefficient is of O {9em^/MQ CD ) and is suppressed by two powers of tu^/Mqcd with respect 
to go- The second term in Eq. (|2"4"]l is basically a correction to go of the same order. Of the TV operators 
with two derivatives, fj — O (Sml/MQ CD J represents a contribution to the radius of the TV pion-nucleon 
form factor, while the remaining operators relevant at this order are relativistic corrections to the go term 
with coefficients fixed by Lorentz invariance. As in the case of the A = 1,2 TV Lagrangians, the TV 
coefficients in Eq. (|24p are related to isospin-violating parameters in the A = 3 TC Lagrangian, which 
are at present poorly determined (see discussion in Ref. [25]). 

For the nuclear potential we will need in addition operators involving more nucleon fields, which are 
constructed in the same way as were the / = 0,2 interactions in Ref. [25 . In the / = 4 sector of the 
theory, the first contribution to the Lagrangian comes at A = 2, 

4=4,74 = —jT])* ' (^NtN nn + a^NtS^n ns"n) , (25) 

in terms of two TV parameters 7 S)0 -. Just as for / < 2, here too there is a link with isospin-violating 
operators, in this case 



(2) _ 1> AT L. _ * . A AT AT AT 0„ AT _ *L 



4=4,t, = fN [T3 - j^n -t)NNN + 2 la N \t 3 - -rjS^N NS»N, (26) 

which generate the dominant contributions to the short-range isospin-violating two-nucleon potential 
[64l [51]. The isospin-violating coefficients 7» = 0{ern^l 'F%Mq CD ) can be seen as low-energy remnants 
of p-ui mixing [HI] and ax- fx mixing [77]. The TV parameters are related to them by 7; = pji = 
O(0rn^/ F^Mq CD ). The operators in Eq. ([2"5"j) are not relevant for the nuclear potential to the order we 
work, but contribute to the three-nucleon TV potential at next order. 

The first operators relevant to the calculation of the TV two-nucleon potential are 

4=4,74 = ( X - J575) [CiNNd^NS-N) + C 2 NtN ■ V^NS^tN)] , (27) 
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where the Cj are two new TV parameters. These interactions are related to isospin-breaking operators 



4=4,t/ = [CxNNd^NS'N) + C 2 NtN ■ V^NS^fN)} , (28) 

with coefficients d = 0(eml/F%M^ CD ). The TV parameters are C, = pd = 0(9ml/F%M^ CD ) 
and, therefore, contribute to the TV potential at order v — 3. The coefficients Cj could in principle be 
determined from pion production in the two-nucleon system and/or from isospin- violating three- nucleon 
forces. However, even lower-order isospin- violating three-nucleon forces are very small [55] . so prospects 
for extracting Ci from TC data are grim. As before, we do not write in Eq. ([27]) other interactions that 
contain more pion fields. 

One can continue the construction of the TV Lagrangian not only to higher orders, but also to more 
nucleon fields. For any given /, the dominant TV terms are expected to be those without derivatives that 
transform like the fourth component P4 , just as in Eq. (|25[) . This type of term will involve an odd number 
of pions in addition to the / nucleon fields. It thus contributes at tree level only to (/ + 2) /2-nucleon 
forces, or to absorption/production/scattering of pions on //2-nucleon systems. The first short-range 
//2-nucleon TV force comes from four- vectors P that involve one derivative, as in Eq. (|2"T|) . Since nuclear 
forces tend to become less important as / increases (see, for example, Ref. [H]), ^ is unlikely that terms 
with / > 6 need to be constructed. The first / = 6 operators that contribute to the three-nucleon TV 
potential appear at next chiral order, A = 4. 



2.3 T- Violating Chiral Lagrangian Prom Dimension-6 Sources 

The smallness of 9 leaves room for other sources of T violation in the strong interactions, which have 
their origin in an ultraviolet-complete theory at a high-energy scale, such as, for example, supersymmetric 
extensions of the Standard Model [TT] . Well below the scale Mp characteristic of T violation, we expect 
TV effects to be captured by the lowest-dimension interactions among Standard Model fields that respect 
the theory's SU c (i) x SUl(2) x fTV(l) gauge symmetry. Just above Mqcd, strong interactions are 
described by the most general Lagrangian with Lorentz, and color and electromagnetic gauge invariance 
among the lightest quarks, gluons, and photons. The effectively dimension-6 TV terms at this scale can 
be written as [8j [9j [TOl 03] 

C T6 = - l -q (d + d 3 r 3 ) o^ifq - lq (d + d 3 r 3 ) a^ 5 X a q G% 
+ ^ £ ^ f a b c G a pG b PGla + l TmEi (qqgrfg _ 0q . ^5 ? ) 

+ ilm£ 8 (q\ a q q\ a rfq - q\ a fq ■ q\ a fz 7 5 9 ) , (29) 

in terms of the photon and gluon field strengths F^ and G" , the standard products of gamma matrices 
7 5 and a^ v in spin space, the totally antisymmetric symbol e^ x<T , the Pauli matrix Tj in isospin space, 
the Gell-Mann matrices A a in color space, and the associated Gell-Mann coefficients f abc . 

In Eq. (|2T)|) the first (second) term represents the isoscalar do (do) and isovector ^3 (d 3 ) components 
of the qEDM (qCEDM). Although these interactions have canonical dimension 5, they originate just 
above the Standard Model scale My/ from dimension-6 operators [8] involving in addition the carrier of 
electroweak symmetry breaking, the Higgs field. They are thus proportional to the vacuum expectation 
value of the Higgs field, which we can trade for the ratio of the quark mass to Yukawa coupling, m q / f q . 
Writing the proportionality constant as eS q f g /Ml (4ir5 q f q /Mp), 

do,3 - O (^o,3^0 , rfo,3 ~ O (^0,3^0 , (30) 
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in terms of the average light-quark mass m and the dimensionless factors (5o,3 and (5o,3 representing typical 
values of 5 q and S q . The third term in Eq. ([29)1 [9] is the gCEDM, with coefficient 

in terms of a dimensionless parameter w. The fourth and fifth operators [TU1 [H] are TV FQ operators, 
with coefficients 

in terms of further dimensionless parameters cri.s- The sizes of <5o,3, <^o,3i w, and cr^g depend on the exact 
mechanisms of electroweak and T breaking and on the running to the low energies where non-perturbative 
QCD effects take over. The minimal assumption is that they are 0(1), 0(g s /Air), 0((<7 s /47r) 3 ), and 0(1), 
respectively, with g s the strong coupling constant. However they can be much smaller (when parameters 
encoding TV beyond the Standard Model are small) or much larger (since f q is unnaturally small). For 
discussion and examples, see Refs. [4"I[TT]. 

The dimension-6 operators in Eq. (|29j) have different transformation properties under chiral symmetry. 
The isoscalar and isovector qEDM (qCEDM) transform, respectively, as the fourth and third components 
of two 50(4) vectors V and W (V and W), with 

^ x ( qa^fq \ p TJ/ _ 1 ( -iqa^j 5 rq 

and 

~ 2 ^ iqa^-f 5 X a q J^^' = 2 ^ qa^X a q 

In contrast, the gCEDM and the two TV FQ operators S^g are singlets of the chiral group, respectively 

iw = \e^f abc Gl p G b /Gl a (35) 



V = l A i^q)^ W -"{ >i ^' (33) 



-. 1 , qa^rX a q \ w = lf -iQ^7 5 r\ a q . 



and 



J ggl = ~ (gg qi^q - qfq ■ qrij 5 q) , (36) 

I qqS = - A (q\ a qq\ a l^q~q\ a Tq-q\ a Tl^q). (37) 

They do not break chiral symmetry. 

The different chiral properties of various TV sources have profound implications for the form and rela- 
tive importance of nucleon-pion and nucleon-nucleon TV couplings in the effective Lagrangian. Effective 
interactions are constructed to transform in the same way as the sources at quark/gluon level. Thus, W3 
leads to interactions proportional to £3 that only appear from in the tensor product between P4 and 
P3 from the quark-mass-difference term, and are thus proportional to em 2 9 (for example, the g~\ term in 
Eq. 421). 

On the other hand, V4 generates exactly the same interactions as P4, so So and 9 contribute similarly 
to low-energy observables. For qEDM, hadronic interactions arise from integrating out at least one hard 
photon, which leads to further breaking of chiral symmetry in the form of tensor products of V and W 
with a antisymmetric chiral tensor |64j . The contributions of the qEDM to purely hadronic couplings, 
like pion-nucleon or nucleon-nucleon couplings, are suppressed by the electromagnetic coupling constant, 
a cm /47r ~ em^ r /(2TTF 7r ) 3 . electromagnetic currents. We do not explicitly construct the TV potential 
from the qEDM. Because they are chiral invariant, the gCEDM and the two TV FQ operators lead to 
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exactly the same effective interactions, although, of course, with different strengths. The interactions 
from gCEDM are the hadronic matrix elements of Iy/ and are proportional to w in Eq. (l3~Tj) , while 
those from the TV FQ operators are the hadronic matrix elements of I qq i in Eqs. (|36|) and (l37l) and are 
proportional to <7i in Eq. (|32[) . For simplicity of notation, in the following w denotes both w and <ii%, 



{ty, (71,(78} w. 



(38) 



For none of the dimension-6 sources there is an immediate, useful connection to TC operators as for 6. 

The pion-nucleon Lagrangian from dimension-6 sources was constructed in detail in Ref. |26j . The 
only terms we need in the following belong to the lowest-order Lagrangian, 



L 'f=2,?6 



In this Lagrangian: 



go 
F„D 



n ■ NtN - 



9i 
F„D 



7T 3 NN ■ 



!I2 



F 2 



(v 



Djt x D^n) ■ NS^tN - 



T3 



V\ ■ DlTT 



FID 
NtN. 



(T7 3 TT • T 7T r 3 



N 



(39) 



For qCEDM, A„ = — 1. At this order the qCEDM contributes only to the non-derivative couplings 
go and g\. In contrast to the 9 case (|24|) . the elimination of the pion tadpoles induced by the 
isovector component W 3 of the qCEDM generates an interaction of exactly the same chiral form as 
other existing interactions, and for simplicity we absorb the tadpole contribution in the gi. Then 



go = O{{5 + sS 3 )mlM QC D/Mp and g x 



0(S 3 mlM QCD /M 2 ) 



• For qEDM, A T = 2. At this order the qEDM contributes also to the non-derivative coupling 
§2, which arises from the tensor product of W 3 with the antisymmetric chiral tensor generated 
by a hard photon. In this case, again after tadpole extermination, go = 0(a e m(c>o + ^3)(1 
,2 m__ ,a„ia*\ ^ - 0{a cm {So + 6 3 )mlM QC D/4TrM 2 ), and g 2 = 0{c 



e)mlM Q c D /^M^), g x 



[a cm 6 3 mlM QC D/^M^). 



For chiral-invariant sources, A T = — 1. It is not possible to construct chiral-invariant, TV pion- 
nucleon couplings with zero or one derivatives. The leading TV couplings then contain either 
two derivatives or one insertion of the quark mass, and the couplings <?o,i and «o,i all appear at 
the same order. More precisely, after the elimination of the pion tadpole, we find go — 0(w(l + 
e 2 )mlM QC D/M 2 ) and g x = 0(wemlM QC D I 'M%), while z ,i = 0{wM QCD /M 2 ). 



In the case of qEDM, we need also the photon-nucleon TV interactions 



■7=2,^6 



2N 



d 1 



2n 2 \ 

fJd) 



di 



T3 



2tt 3 . 
F 2 D' 



S^Nv'F^, 



with do,i = 0{e8o :3 m 2 T / M^,Mqc d) short-range contributions to the nucleon EDM. 

The various dimension-6 sources lead to different interactions also in the / > 4 sectors, 
chiral-invariant sources, there are further leading-order interactions with / = 4, 



C 



(-1) 



dNNd^NS^N) + C 2 NtN ■ V^NtS^N), 



(40) 



For the 



(41) 



with coefficients C1.2 = 0(wMqcd/F 2 M^,). Note that because of different chiral-transformation prop- 
erties, the form of Eq. (l4"Tj) is different from the contact interactions stemming from 9, Eq. (f2~T|) . For the 
remaining dimension-6 sources, the form of contact interactions can be different still. They appear from 
tensor products of symmetry-breaking operators and are subleading. The qCEDM and qEDM generate 
operators at chiral index A = A T + 1, which have schematically the same form as Eq. (|25|) . but a richer 
isospin structure. These operators do not contribute to the two-nucleon potential at tree level, but they 
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are relevant for three-nucleon forces just beyond the order we consider. Contributions from these sources 
to the two-nucleon potential are found at A = A^ + 2, 

4t£j$ = ( X - U) l°^ N 0*(XS»N) + C 2 NtN ■ V^NS^tN)} 

•2 



1 - [CuffnNdp (NS^N) + C 6 NNV^ (N n S»N)] 



(<*3k ~ Jf^*) Nr k NV^ (NS^nN) 



2 

*13 + (7T37TJ -7r 2 <5 3/ ) 



(42) 



where the Ci are new coefficients. The qCEDM only contributes to Ci.2,3.4 in lowest order, with C1.2 = 
C((<5 + eS 3 )ml/F^M^M QC D) and £3,4 = 0{5 3 ral/ F^M^M QC d)- The qEDM contributes also to 
^5,6,71 because of tensor products with the antisymmetric tensor generated by a hard photon, and we 
have C x , 2 = 0(a cm (S Q + 5 3 )(1 + e)m 2 jAnF^M^M QC D), C 3A = O(a em (S + Ss^/AttF^M^Mqcd), 
C 5 , 6 = O{a cm S ml/ATTF^Mf,M QC D), and C 7 = 0(a cm S 3 m^/ AttF^M^Mqcd))- Comparing Eq. (02J 
to Eq. (|27p. we once again see that, differently from the 9 term, the qCEDM and qEDM generate TV 
and isospin-breaking operators of the same importance as isoscalar TV operators. The four-nucleon 
operators in Eq. (l4"2l have, this time just like the 6 term, chiral index that is two units bigger than that 
of the leading pion-nucleon Lagrangian. As a consequence, short-range contributions to the two-nucleon 
potential from these sources arise only at next-to-next-to-leading order. Such high orders will not be 
considered explicitly below. 



3 Pionless Theory 

Before we discuss the TV potential in ChPT, it is instructive to consider a much simpler EFT. At 
momenta much smaller than the pion mass, pion degrees of freedom can be integrated out and one is left 
with a pionless EFT, in which the interactions are represented by operators involving only nucleon fields. 
If we denote by M nuc ~ 100 MeV the scales associated with pion physics, this EFT applies to processes 
where all momenta Q <C M nuc . Power counting in this EFT is different from ChPT and is reviewed in 
Ref. 0SJ. 

The lowest-order TC two-nucleon interactions can be taken as [66] 

c c 

Cf' T = ^-NNNN ^-NtN ■ NtN, (43) 

where Cn jTT = 0(47r/mArH), with H < M nuc a low-energy scale. The corresponding potential in momen- 
tum space is simply 



(44) 



where r^/2 is the isospin of nucleon i. These interactions affect only the two S waves. Since the effect 
of free two-nucleon propagation is ~ mpfQ/Air, for momenta Q > H these interactions have to be iterated 
to all orders [TSj [7J5] ■ Using dimensional regularization with power divergence subtraction !79 : at a scale 

4tt / 1 \ _1 4tt / 1 \ _1 

Cos = C n - 3C TT = — -M , C Qt = C u + C TT = , (45) 

m N \a s J m N \a t ) 

in terms of the isospin-singlet ( 3 Si) and -triplet ( 1 <So) scattering lengths, a s and a t . Because the coeffi- 
cients Cn. rr subsume physics at the scale of the pion mass, their scaling is different from the one in the 
pionful EFT. 
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In leading order, the 6 term and the dimension-6 sources induce TV four-nucleon operators similar to 
those in Eqs. (1271). (HT j) and 

£f,r = C U NN dn(NS"N) + CrrNTN -d^NS^rN) 

+C 31 Nt 3 N d^NS^N) + C 13 NN ■ d^NS^rsN) + O^N^Nd^NS^N), (46) 

where Cii )T T,X3,3i,33 are new short-range parameters. In momentum space the interaction Hamiltonian 
is given by 

-|(C 13 -C 31 ) (r 3 (1) -rf) + <?<'>) (47) 



where <? ^ /2 is the spin of nucleon z and q = pi — p[ = p' 2 ~ P2 is the momentum transfer. 

In nucleon-nucleon scattering, operators that break P and T induce mixing between waves of different 
parity. At low energy, the most relevant effect is the mixing between S and P waves, and indeed the 
single momentum in Eq. (1471) can only connect an S to a P wave. At leading order, the P wave is free. 
Since the short-range TV potential involves one S wave, we expect [IS] that in the pionless EFT the 
coefficients CV/ scale as Indeed, the amplitude for a nucleon-nucleon transition can be computed 

from Eq. ((47]) as done in the PV case in Refs. [56, 60 . In leading order, it involves one insertion of the 
TV operators Cy, dressed by the all-order iteration of the appropriate S-wave operator, Cos or Cot. The 
renormalization-group invariance of the amplitude implies that the dj follow a renormalization-group 
equation of the form d(Cij/Co)/dliifi — 0, which is satisfied if the five independent parameters are taken 
to be 

C\ s = On — 3C TT = — — ( fi J , On — On + TT = — — ( fj, 

mjv \a s J m N \a t 

r -f f - 4 ^ 3s ( 1 „V 1 r -r +r 4 ™ 3t ( 1 

<^3s — t>13 — ^31 — AM > ^3* — W3 + <^31 



m a? \a s / mpf \at 

S3 = ^(±-»y\ (48) 
tojv \a t J 

in terms of five /i-independent coefficients c Sj t, 3s , 3 t, 33 t- As in the TC sector, the scaling of the short-range 
parameters is different in the pionless EFT than in ChPT. We can write 

Cu,rr = o(-^c a ,t), Oi 3 , 3 i = o(^-c 3s , 3t ), O 33 =o(^-c 33t ). (49) 
\m N s< J \toatH / \tojvH / 

In order to estimate the coefficients c S:t ,3 S ,3 tj 33 t , we use naive dimensional analysis |17[ 1651 |S] with 
the pionful EFT as the underlying theory. We then find that at leading order the isoscalar c s t receive 
contributions from all the sources, 



.,-> i , ,- . 5 sMqcd cn cm ,Mq C d Mqcd . 

while the isospin-breaking c 3s ^t only from the dimension-6 sources 



m I x m Qcd "cm,, . x \ Mqcd Mqcd \ , K1 , 

C3s < 3t = ° (^isp ^r {So + ^^p^^r ) ' (51) 
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and C33 t only from the qEDM, 



c 33f - O 



o 



.■m c MQCD 

17 3 "mF" 



(52) 



In general, one would expect five possible amplitudes connecting S to P waves [3H US] : three — one 
for each possible value of 73 = 1,0,-1 — to describe the mixing of the isotriplet 1 Sq and 3 Po waves, 
one for the mixing of the isosinglet 3 Si and 1 P\ states, and one for the mixing of nucleons in the 3 Si 
configuration with the isotriplet 3 Pi wave. The 9 term yields a short-range potential in the form of the 
isospin-conserving terms of Ref. [34]. Because the TV operator in Eq. (|21|) is isoscalar and isospin 
violation is a subleading effect in ChPT, for which the pionless EFT is the low-energy limit, the 9 term 
does not contribute at leading order to quantities that violate both T and isospin. The two terms 
contribute to 3 S\- l P\ mixing and to 1 S$- 3 Pq mixing, in equal way for the three I3 configurations. The 
3 S\- 3 P\ mixing vanishes at leading order, a fact that has important consequences for the estimate of the 
deuteron EDM 28 . If S3 and e are different from zero, the qCEDM and the chiral-invariant TV sources 
also contribute to isospin-breaking TV observables at leading order. The operator Czt is proportional to 
the third component of the total isospin of the two-nucleon pair, and thus it does contribute to 1 Sq- 3 Pq 
mixing, but only for I3 = ±1. Cj, s is instead proportional to the total spin of the two nucleons, and it 
is relevant to 3 S\- 3 P\ mixing, and, consequently, to the deuteron EDM. Only the qEDM produces full 
isospin breaking in leading order. 

A potential with just these five short-range terms was considered recently |47j . Five low-energy quan- 
tities — such as the spin rotation of a polarized beam and the longitudinal polarization of an unpolarized 
incident beam in neutron scattering on a proton target [46] at different energies — are needed to determine 
the parameters c a ,t, c 3s ^ t , and c 33t . 

There is, however, an important extra ingredient that needs to be added at low energies: one-photon 
exchange where one of the vertices originates in the nucleon EDM, as in Eq. ([40]) . see Fig. [1] This long- 
range potential is particularly important for qEDM, since for this source the short-range interactions 
have suppression by a em /47r from the hard photon, as can be seen in Eqs. (T5U|) . (p)T]) . and (|5"2"|) . We find 



ie 

T 

ie 



(do + d,) (r 3 (1) 



>> 






A 2) ) 




-a( 2 



+ (di - d ) (t. 



-(2) 



(53) 



where 



do,i = 



M^Mqcd 



(54) 



are the isoscalar and isovector components of the nucleon EDM. In principle, the coefficients do,i in the 
pionless theory are different from the short-range contributions to the nucleon EDM in the pionful theory 
introduced in Eq. (14"0"|) . However, since in the case of TV from the qEDM pion-loop contributions to the 
nucleon EDM are suppressed [51], at leading order they exactly match. 

The photon-exchange potential has the same spin/isospin components of the short-range potential 
(|47|) . but it acts on all partial waves. While for the other sources this type of potential is important only 
for very low momenta, for qEDM it is enhanced with respect to the short-range potential (47) by a factor 
of M ntic /H for Q ~ K, and it dominates throughout the regime of the pionless EFT. 

In subleading orders more derivatives appear in two-nucleon contact interactions and photon exchange. 
An important issue is the order where few-nucleon forces first appear. On the basis of naive dimensional 
analysis, we expect them to be also of subleading order. In the PV TC case, this is confirmed for the 
three- nucleon force by a more detailed analysis based on renormalization-group invariance |58| . Since the 
powers of momenta involved here are the same, the same conclusion should hold. The potentials (|47p 
and (1551) should then be sufficient for most TV applications of the pionless EFT. 
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Figure 1: One-photon-exchange diagram contributing to the long-range TV two-nucleon potential. The 
solid and wavy lines represent nucleon and photon, respectively; a square stands for the TV photon- 
nucleon couplings in C^J_ 2 p 6 (|40|) , while the filled circle represents an interaction from £y< 2 t C3)- Only 
one possible ordering is shown. 

4 The TV Potential in Momentum Space 

In processes involving momenta Q ~ M nuc , which presumably comprise the bound states of most nuclei, 
pion effects are important and pion degrees of freedom should be included explicitly in the theory. In 
this section we use the interactions given in Sect. [2] to compute the TV nuclear potential in momentum 
space for the TV sources of dimension up to 6. 

In the lowest orders, the TV nuclear potential involves only two nucleons. We write the two incoming 
momenta as p\ — P/2 + p and p2 = P/2 — p, and the two outgoing momenta as p[ = P/2 + p' and 
P2 = P/2 — p'. The TV potential in momentum space can be expressed as function not only of the 
momentum transfer q = p — p', but also of the center-of-mass (CM) momentum P and of the variable 
K = (p+p')/2: Vf = Vp(q, K, P). Expressions for the potential in the CM frame are obtained by setting 
P = 0. Notice that although some of the terms below vanish in the CM frame, they can be relevant to 
the calculation of the TV electromagnetic form factors of deuteron, or for calculations of T violation in 
nuclei with A > 2, where the interaction with the photon or other nucleons changes the CM momentum 
of the nucleon pair. 

4.1 9 Term 

In leading order, the 0-term nuclear potential comes from the OPE diagrams of Fig. [5J with TC and TV 
pion- nucleon interactions taken from T and £^jl 2 fi m Eqs. (|15j) and (|22[) . respectively. The strong- 
interaction vertex introduces a factor of gAQ/F-K, while the TV vertex brings in a factor go oc iti^/Mqcd- 
As a result this contribution goes as Mq CD and it is of order v = 1. In momentum space, the expression 
for the potential is simply 

yd) = ^(i) . f (2) (V(D _ . _JL_, (55) 

which agrees with Ref. [29] . Just like the 9 term contribution to the potential (j47|) in the pionless theory, 
this OPE potential contributes to 1 Sq- 3 Pq and 3 S'i- 1 Pi mixing, but not to the isospin-violating 3 S\- 3 Pi 
mixing. At this order, there is a single unknown TV parameter, go . Contrary to the PC, TC case [JS] 
and more like the PV, TC potential [56], pion physics is enhanced relative to short-range physics due to 
the absence of a derivative in the simplest pion-nucleon TV interaction and the presence of one in the 
simplest TV two-nucleon contact interaction. 

Since in nuclei the OPE from the 1 = pion-nucleon coupling go is suppressed by the factor (N — Z) /A, 
it is interesting to pursue higher orders, up until the 1=1 pion-nucleon coupling gi, whose OPE is 
not affected by such suppression, appears. This means up to v = 3, that is, including corrections of 
0(Q 2 /Mq CD ) with respect to the leading TV potential. According to Eq. (TT3l) . corrections at orders 

v = 2, 3 come from one-loop diagrams involving £y°< 2 t an< ^ ^7=2 fi only, and from tree diagrams with 
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Figure 2: OPE diagram contributing to the leading TV two-nucleon potential. The solid and dashed 
lines represent nucleon and pion, respectively; a square stands for the TV pion-nucleon couplings go in 
Cjl 2 fi or 5o,i,2 and %\ in Cj^vTe P^- wnne the filled circle represents an interaction from £y< 2 t 
(fT5|) . Only one possible ordering is shown. 
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Figure 3: Box, crossed, and triangle TPE diagrams contributing to the subleading TV two-nucleon 
potential. Notation as in Fig. [2] Only one possible ordering per topology is shown. 



insertions of higher-order terms. The tree contributions come from the four-nucleon TV operators in 
£fl 4 p 4 , Eq. (f27|) . and from OPE diagrams in which either the TC or the TV vertices originate in the 
power-suppressed / < 2 Lagrangians. 

The most important loop diagrams are from TPE, depicted in Fig. [3J The T-odd pion-nucleon coupling 
go and one of the strong-interaction vertices bring in a factor of gogA/F 2 . The other two vertices of the box 
and crossed diagrams of Fig. [3J are strong-interaction pion-nucleon vertices from Eq. (|T5j) , and combined 
with the (47r) 2 from the loop integration, they yield the suppression factor g\j (AttF^) 2 ~ 1/Mq CD . For 
the triangle diagrams, the seagull vertex is the Weinberg- Tomozawa term also from Eq. (|15[) . which 
brings in a factor of 1/F% that, combined with the (47r) 2 from the loop, also leads to a suppression of 
l/(47rF 7r ) 2 ~ 1/Mq CD . All these diagrams are thus of order Mqq D . Care of course has to be taken with 
the subtraction from the box diagrams in Fig. [3] of the iterated static OPE, which is infrared enhanced 
and already included in the computation of wave functions at lower order. Following the procedure 
described for instance in Ref. [S3], the subtraction is accomplished by exploiting the identity 

2n5(v ■ k). (56) 



— v ■ k + is v ■ k + is 

When Eq. ([56]) is used in place of one of the nucleon propagators in the box diagrams, the first term 
on the right-hand side leads to a contour integral over the Oth component of the loop momentum, which 
can be performed without picking up the nucleon poles and is free of the infrared enhancement discussed 
in Sect. [5J while the delta function corresponds to the two-nucleon pole and must be discarded in the 
calculation of the potential. For the crossed-box and triangle diagrams, instead, it is always possible to 
avoid the nucleon poles, and these diagrams only contribute to the potential. 

The TPE diagrams in Fig. [3Jare ultraviolet divergent. We regulate them in dimensional regularization 



16 



in d spacetime dimensions, where divergences get encoded in the factor 



In 4-7T, 



(57) 



where je is the Euler constant. We denote by [i the renormalization scale. Proper renormalization 
requires that sufficiently many counterterms appear at the same order to compensate for the L and /i 
dependence of the loops. Indeed, here this dependence can be absorbed by the renormalization of the 
contact interaction C 2 from £y= 4 y 4 , Eq. (j2"7j) , which we do by redefining it through 



Co Co 



1) [L + In ■ 



2 (ffi - 1) 



(58) 



Note that we chose to absorb in C 2 some finite constant pieces. TPE does not renormalize the coupling 
Ci at this order. With this redefinition, the contact interactions yield the short-range potential 



V- {3) (q) = -- 



C 1 + c 2 f^-r^] (jO-jP)).,. 



(59) 



which is formally identical to the leading 8 potential in the pionless EFT, Eq. (|4"7|). The couplings, 
however, are different. We can see from Eq. (|58|) that the natural size of the coefficients Cj is, as 
advertised, #m 2 / 'F^Mq CD , implying a suppression of Q 2 /Mq CD with respect to TV OPE. 

Once the divergent, short-range part of TPE has been lumped with the contact terms, we are left with 
the non-analytic contributions of medium range, 



(3) 



. ZgogA 
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V 4m 2 
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(60) 



in terms of the functions 



T(x) 



■ In (Vi + VT+Ir) 



-B(a:). 



(61) 



As the leading OPE potential, Eq. (l55t , the TPE potential is a function only of the momentum transfer 
q. The scale of momentum variation is, as one would expect, 2777,^. TPE and leading OPE share the same 
spin-isospin structure, which means they can only be separated if we probe their different momentum 
dependences. 

A much richer structure arises from the remaining v < 3 contributions to the two-nucleon TV potential, 
which come from the OPE diagrams depicted in Fig. |4j Doubly-circled vertices in the first two diagrams 
denote 0{Q 2 /M^ CD ) corrections to the TC and TV pion-nucleon couplings, given by the operators in the 

Lagrangians £y< 2 ^ T /> ^/<2T/' an< ^ ^7=2 ;T4 f° un d in Eqs. (fTB]) . (fT7|) . ([2"D]). and (j2"4"|) . The last diagram is 

— . . (2) (1 2) 

proportional to corrections to the pion mass in £j< 2 TI and £^< 2 and to the nucleon mass difference 

(12) _ -1 

in £y-< 2 7y Note that, as we argue shortly, there are no further loop diagrams to consider explicitly. 

Corrections that originate in the pion mass are closely connected to the leading OPE, Eq. (|55l) . At the 
order we are considering, the pion mass receives corrections from one-loop diagrams, which we absorb [18] 
in the renormalization of the coupling Am 2 in Eq. (fTTl) . With the definitions of Eqs. dl~5|) . (fTT|) . (ITS)) , and 
(Til?]) , the physical masses of the neutral and charged pions are, respectively, m 2 = + + <5m 2 — 
(135 MeV) 2 and m 2 ± = m 2 + Am 2 . + <5m 2 = (139.6 MeV) 2 [72]. The isospin-symmetric correction to 
the pion mass can be accounted at v = 3 by substituting m 2 — > m 2 + Am 2 in the leading-order TV 
potential. Isospin-breaking corrections come from the different masses of the neutral and charged pions. 
With the assumption a cm /47r ~ £to^/Mq C£) , which is numerically reasonable, the pion mass splitting is 

dominated by the electromagnetic contribution (5m 2 , which gives rise to a potential of order v = 2. The 
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Figure 4: OPE corrections to the #-term two-nucleon potential up to order 0(Q 3 /Mq CD ). The double 
circles denote vertices in the A = 1,2 TC chiral Lagrangians, £j< 2 ti P^l - £/<2 tj C3> anc ^ ^/<2 T/ 
(|20|) . The doubly-circled square denotes vertices from the A = 3 TV Lagrangian, 

£fl 2 ,r4 ®- other 

notation as in Fig. [51 Only one possible ordering is shown. 



quark-mass-difference contribution Sm^ appears at v — 3, when one should also consider diagrams with 
two insertions of 5m„. The sum of these components generates two structures, an isoscalar 

y t ^ . _^ ( 2 J m > + , mi _ 2 jM|L) .... . ? « (,<» _ . _J_ m 
and an isotensor 

Cm - 'If (** - *»i - - ?w - ?<2) ) (*" - s,2) ) ■ (Fw < 63 > 

These isoscalar and isotensor components can be rewritten using the physical pion masses. The isoscalar 
component, Eq. (|62[) . can be obtained by using the physical values of the neutral and charged pion mass 
in the leading potential, Eq. (|55j) . We can write 

V^\q) + V^\q) = igj ( ^ 2 2 + ^ 2 ) f CD • r W (V« - a' 2 >) ■ g, (64) 
9 ^ ' 9 < a ^ ' 3F 2 \q2+ml ± q 2 +m 2 ) \ J 

which, expanding in Sm^ and Sm 2 , reproduces Eq. (|62p . The combination of neutral and charged pion 
propagators in Eq. (|64p represents an "average" pion static propagator, which naturally appears in the 
isoscalar contribution. Similarly, we can rewrite the tensor component as 

vir-w - -<!§ - (*w - ■ **) (*» - *») ■ * w 

In applications to nucleon-nucleon scattering, this tensor component would contribute at low energies to 
1 Sq- 3 Pq mixing, affecting proton-proton and neutron-neutron (I3 = ±1) and neutron-proton (^3 = 0) 
scattering differently. It is worth stressing that, in contrast to phenomenological approaches, the isotensor 
component coming from the 9 term is not a leading contribution. 

Corrections from the nucleon mass come in several guises. The use of a heavy-nucleon field ensures 
that the large scale appears always in denominators. In the isospin-symmetric limit the first effects 
of m N enter in the A = 1. 2 TC Lagrangians, Eqs. (JTSJl and fI7j). the A = 3 TV Lagrangian, Eq. 
and, via the on-shell condition for the nucleons, the energy of the potential-pion propagator. They yield 
relativistic corrections to the leading OPE with the same spin-isospin structure, 

V^(q, K, P) = (k 2 + ^- 7^41 r W • r ^ (>) - a^) ■ (66) 

e ' c y ' F 2 m 2 N \ 4 4 q 2 + m 2 J V / q 2 + m 2 
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and with new structures, 



V^(q,K,P) = 



= , 9A9o f(i) f(2) 1 

^F 2 m 2 N ' q 2 + m 2 



P-q 



(V(i)_5K2)).| + (>) +a K2)) 





(67) 



Note that in writing Eqs. (|66l) and (|67p we have omitted pieces that vanish due to energy-momentum 
conservation for on-shell nucleons, which implies K ■ q = p 2 — p' 2 = 0. This potential includes the 
contribution of the 1/mjy correction to gA in Eq. (|16|) . Naively, one would expect this correction to 
contribute at order v = 2; however the interaction brings in a factor of v • q, which, for on-shell nucleons, 
becomes v ■ q = q - P/2m N , suppressing the potential by a further factor of l/m N . There is a subtlety in 
this argument. When one performs the integral involving both OPE with a pion energy in the numerator 
and another interaction in the potential, and picks the pion pole, one gets a one-loop contribution to the 
potential. However, by power counting, such diagrams are suppressed by a further Q/Mqcd- (This of 
course does not preclude enhancements by factors of it that in principle might affect any ChPT loop, but 
are hard to incorporate in power counting.) 

Corrections to the nucleon mass can be removed from nucleon propagators by redefinitions of the 
nucleon field. The chiral-symmetry-breaking correction to the nucleon mass, Amjv in Eq. (1161) . can be 
absorbed in rnjv, m N — >■ rnjv — Arnjv, by a redefinition of the nucleon field of the same type of that which 
eliminates the mass from Eq. (p~5|) in the first place. The isospin- violating nucleon mass splittings 5m^ 
and 8mN can be dealt with the field redefinition of Ref. [54], which leads to Eq. (|20|). The corresponding 
potential linear in 5m n has a 1/m^ factor, 



V^(q,K,P) 



9o9A Sm N 



f « x r ( 2 ) 



F 2 m N q 2 + m 2 . 



q 2 + m 2 



P-q 



This potential has the right quantum numbers to produce 3 £i- 3 Pi mixing, and therefore must be included 
in a calculation of the deuteron EDM. In addition, there are terms quadratic in Sm^, which generate an 
additional contribution to the isoscalar and tensor potentials in Eqs. (f6"2"]) and (f6"5)) . 



■ 9ogA x 2 
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(69) 



Note that the electromagnetic correction to the nucleon mass, Sm^, does not appear in the equations 
above. The reason is that in Eq. (|20p it appears with a v ■ drf factor, which again brings a v ■ q and 
consequently an extra suppression ~ Q/tojv- The contribution from Sm^ is thus next order. Up to such 
higher-order terms, we can make Sm^ — > SmN + Stun = m n — m p in the expressions above. 

Finally we arrive at contributions from v — 3 effects in the pion-nucleon vertices. At this order, several 
contributions can be absorbed into redefinitions of the couplings gA, go, and Ci. One-loop corrections 
to gA do not introduce any non-analytic contribution and, for an on-shell nucleon, they renormalize the 
coupling d,A in Eq. (TIT]) [80 . The operator with coefficient <1a gives rise to a potential like Eq. ([55]) . with 
gA replaced by — gA^A- For simplicity we absorb <1a in gA, gA ^ .9a(1 + Similarly, the calculation of 
the pion-nucleon TV form factor in Ref. [3S] shows that the one-loop corrections to g do not introduce 
any non-trivial momentum dependence, so they simply renormalize the coupling A^ in Eq. (I24[) . These 
ml corrections to g can be absorbed in it, go — > go — Ag + g^dml/ml. As for the operators with 
coefficients ca in Eq. (|17l) and fj in Eq. they give potentials of the form 



f(l).f(2) 



?(2) 



_^U).f(2) 



f(2) 
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■q, (70) 
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which cannot be distinguished from those of <M3o and C_, and can therefore be absorbed into further 
redefinitions of Qa, go, and C_. Now the Goldberger-Treiman relation for the strong pion-nucleon constant, 
9-nNN — Zm-NgA/ 'Fn, applies without an explicit correction. If for the pion-nucleon coupling constant we 
use girNN = 13.07 [81], then in the leading-order TV potential we should use gA = 1-29. 

The remaining contributions come from vertex corrections, both in the TV sector via the TV pion- 
nucleon coupling tt^NN in Eq. (|24|). and in the TC sector via the isospin-breaking pion-nucleon axial- 
vector coupling d^NS^N in Eq. (gDJ). We find 

+ (SAS, + S -f) (r„" - rf>) (JO + gm) ■ f} , (71) 

where we redefined g~\ to absorb the tadpole contribution g~\ — > g~\ + 2<7oATOjv<Jm„ / 'Sm^m^. The first 
structure contributes to 1 Sq- 3 Pq mixing. Being proportional to I 3 , the contribution vanishes in the case of 
neutron-proton scattering, and is only relevant for proton-proton or neutron-neutron scattering. Because 
of its isospin structure, it does not affect the 3 Si- 1 Pi and 3 Si - 3 Pi channels, and, in particular, it is not 
relevant for the calculation of the deuteron EDM. The second structure, in contrast, contributes to 3 Si- 
3 Pi mixing, and, consequently, to the deuteron EDM. Its contribution vanishes in the other low-energy 
channels. 

Note that loop diagrams involving the leading S-wave TC four-nucleon operators and a TV pion 
exchange all vanish. The analysis of Refs. [67l [68] showed that some higher-wave TC four-nucleon 
operators are less suppressed than expected on the grounds of naive dimensional analysis, and they must 
be included in the leading-order / = 4 TC Lagrangian. Loop diagrams with P-wave operators and a TV 
pion exchange do not vanish. However, these diagrams do not depend on the momentum transfer q and 
they simply renormalize the couplings C\ and C<z- 

One can proceed in the same manner to construct higher-order potentials. At next order there are 
further OPE and TPE, and also one-photon-exchange, contributions to the two-nucleon potential. There 
is also the appearance of the lowest-order three-nucleon TV potential, which arises from essentially three 
mechanisms: (i) a TPE component cx gAgo/m^F^ involving a pion energy in a Weinberg- Tomozawa 
seagull vertex; (ii) a TPE component cx g A ho/F% involving the seagull vertex from £f}_ 2 fv ^1' P^l) ; 

and (Hi) a one-pion/short-range component cx gAji/F^ involving the short-range pion- two-nucleon inter- 
( 2) _____ 

actions from CfL± 7*4. Eq. (P25l) . The fact that, in the absence of an explicit delta isobar, the three-nucleon 
potential first shows up three orders beyond leading is completely analogous to the TC PC case [50]. An 
important difference is that, because of the relative enhancement of pion exchange compared to short- 
range physics, the leading TV PV three-nucleon force does not include a purely short-range component. 
Thus, this TV PV three-nucleon force is in principle determined by one- and two-nucleon physics. 

4.2 Dimension-6 Sources 

Our attention has been focused so far on the TV potential from the term, in which case the vanishing 
of g~i at leading order makes it important to consider subleading contributions. We now briefly turn our 
attention to the TV potential from dimension-6 sources of T violation. 

For all dimension-6 sources, the leading-order potential contains OPE of the form in Fig. [2] Since 
these sources all generate 1 = and 1=1 pion-nucleon couplings of the same size, we have no motivation 
to go to subleading order in the potential. The leading-order potential, which is a function of the transfer 
momentum q only, should be sufficient for most phenomenological applications. Of course, if needed, 
subleading orders can be derived just as we have done for 9. 

In the case of the qCEDM, the TV couplings go and g\ both appear in the A = — 1 Lagrangian, Eq. 
([39]) . As a consequence, the leading potential from the qCEDM has chiral index v = — 1, and it has 
both an isospin-conserving part, which is identical to Eq. (|55[) . and an isospin-breaking one. The leading 
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potential is 



V, 



+i 



9A9i 
2F2 



.(72) 



For the TV chiral-invariant (CI) sources, that is, gCEDM and TV FQ, the A = -1 TV Lagrangian 
contains both pion-nucleon couplings and four-nucleon operators, see Eqs. ([3"5]l and The addi- 

tional TV pion-nucleon coupling i\ produces a potential of the type ([70]) . and thus can be absorbed in 
redefinitions go go + m~ii and C% — > Ci + 2gA%i/F%. As a consequence, the leading two-nucleon 
potential consists of an isospin-conserving and an isospin-breaking one-pion-exchange contribution and a 
short-distance piece, 
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The qEDM leading-order potential displays further new structures: in addition to the I — 2 pion- 
nucleon coupling §2 in Eq. (1391) . there is also the one-photon-exchange contribution shown in Fig. [TJ 
where one vertex is the short-distance contribution (|4T))) to the nucleon EDM: 
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As for at subleading orders the potential receives corrections from various mechanisms. For all 
dimension-6 sources there are one-loop diagrams with the same topology as in Fig. [31 the square now 
denoting go, g\ or §2- There are also tree-level OPE diagrams of the type in Fig. [H with subleading 
TC and TV pion-nucleon vertices. For qCEDM and qEDM, the one-derivative four-nucleon interactions 
in the Lagrangian (|42]) have to be taken into account, while for TV CI sources one has to include all 
the possible TV four-nucleon operators with three derivatives. For qEDM, one also needs to include 
corrections involving photon exchange. This calculation proceeds along lines that are very similar to 
Sect. O 

In the case of the 9 term, TV three- nucleon forces only appear at NNNLO, one order higher than 
the accuracy of our analysis. For qCEDM and qEDM the situation is similar, but one might wonder 
whether for CI sources, which appear to be more sensitive to short-distance physics, TV three-nucleon 
forces are more relevant. However, also in this case it turns out that the three-nucleon potential is 
a NNNLO effect. The lowest-order three-nucleon potential receives various contributions: (i) a TPE 
component oc gjigo^/mpfF^ involving a pion energy in the Weinberg- Tomozawa vertex; (it) a TPE 
component cx g\%o/mNF^ involving a pion energy in the leading TV seagull %o\ (lit) TPE components 
from TV seagulls in the A = pion-nucleon Lagrangian, which we did not explicitly construct; (iv) a 
one-pion/short-range component cx gjCfijF^, with four-nucleon operators that contain at least one pion 
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field in the A = Lagrangian; and (v) short-range six-nucleon operators, which also appear in the A = 
Lagrangian. From the power counting formula (jl3[) . all these three- nucleon contributions are suppressed 
by three powers oiQ/MQCD with respect to the effects of the two- nucleon potential (|T3[) in the three-body 
system. 



5 The TV Potential in Configuration Space 

The evaluation of T-odd observables in nuclear and atomic systems is often more easily carried out in 
configuration space. In this section we give the TV nuclear potential derived in Sect. |4] in coordinate 
space. We compare this potential with the literature in the next section. 

In the two-body case, it is convenient to introduce the relative position of the two nuclcons r = x\ — X2, 
their CM coordinate X = (x*±+X2)/2, and the conjugate variables — iV r = — id/dr and — iVx = —id/dX. 
Translation invariance constrains the potential to commute with Vx and, therefore, not to depend on 
X , so that in general the potential is a function of r and of the nucleons' relative and CM momenta, 
Vr = Vf(f, V r , Vx)- The relations between the potential in momentum space and in coordinate space 
are detailed in the Appendix. Some care must be taken, and a regularization scheme has to be defined, 
when computing the Fourier transform of functions that blow up as \q\ goes to infinity, as is the case of 
the subleading TV potential. As described in the Appendix, we follow Ref. [82] and define the Fourier 
transform in d dimensions. We apply the ci-dimensional Fourier integration of the momentum-space 
potential before setting d = 4. This method eliminates naturally the divergent factor T(2 — d/2) arising 
from loops and yields a finite result. As we will see, the divergent behavior at large momentum translates 
in a singular ~ \/r potential at short distances. Expressions in configuration space obtained with this 
method are equivalent to the procedure based on old-fashioned perturbation theory |82j . 

In order to write the results of the Fourier transform we introduce a few functions of the magnitude 
r = \ f | of the radial coordinate: 

U(r) = — — [2cxp(-m 7r ±r) + cxp (-m n or)] , (75) 
127rr 

which reduces to the usual Yukawa function U(r) — exp(— m 7r r)/47rr when we ignore the pion mass 
difference; 

W(r) — - — [exp (— ra^±r) — exp (— m„or)l , (76) 
4irr 

which is entirely a consequence of isospin breaking; and the TPE functions 

x[r) = M^k)^r^ (3+3 ^ +/3v)exp( " /3r) ' (7?) 

y W = o (o p ^2 3 [ dx(l+pr) cxp(-/3r), (78) 
with /3 2 = m%/x(l — x). 



5.1 6 Term 

The Fourier transform of the leading OPE potential including the corrections from the pion mass (|64[) 
and from the nucleon kinetic energy (j66]), is 



. (79) 
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where {• • • , • • • } denotes the anticommutator. The remaining pion-mass correction, Eq. (|65p. is 

g| (3r 3 (1) rl 2 ) - ? V • r <»)) - ^) ■ (v r W(r)) , 
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while the Fourier transform of the other relativistic corrections to the leading OPE, Eq. (|67|) . is 
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The nucleon mass-splitting corrections in Eq. become 



while those in Eq. (|69[) read 
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Finally, the Fourier transform of the TPE potential in Eq. (|6T))) reads 

"V, (2 5 2 X(r) - Y(r)) 
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The potential in Eq. (|85[) is singular, and the cutoff dependence it introduces in the evaluation of matrix 
elements and observables is absorbed by the renormalization of Ci in the short-distance potential (I59[) . 
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5.2 Dimension-6 Sources 

At leading order, the potential from dimension-6 sources is local, and depends only on the relative position 
r. 

In the case of the qCEDM, since the potential arises from one-pion exchange, its radial dependence is 
encoded in the Yukawa function U(r), 
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For CI sources, there are additional short-range contributions, 
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The potential from the qEDM is purely long-distance: in addition to pion exchange of range ~ l/m v , 
there is a longer-range component from photon exchange, 
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The various sources thus involve different spin, isospin, and radial dependences. We discuss some 
implications in the next section. 



6 Discussion 

Traditionally the study on T violation in nuclear physics has been carried out by considering the most 
general pion-nucleon non-derivative couplings in a phenomenological TV Lagrangian |31] , which we write 
in our notation as 
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and by inferring from it the TV two-nucleon potential 
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with U(r) defined in Eq. (|75p . When short-distance contributions are included in the model, the most 
general TV two-nucleon local potential with the minimum number of derivatives assumes the form [35] 



,(r) = (a( 1 )-a( 2 )).vL(r)+r( 1 ).r( 2 )V (r) + i(r 3 (1) +r 3 (2) )w 1 (r) 

+ (3ri 1} rf - f W • ?«) V 2 (r)] + i (r,« - rf ) (<x« + a< 2 >) • VV l( r) (92) 



in terms of five radial functions i/o,i(r) an d Vo,i,2(^)- These functions are assumed to originate in one- 
boson exchange [3HG2I1GH]: Pi° n exchange is taken to give long-range contributions to Vo, (Vi +Ui)/2 
and V2, while eta, rho and omega mesons give shorter-range contributions to the same quantities, and to 
Uo and Vi — U\. The five momentum-independent potentials in Eq. (|92[) are treated on the same footing, 
and they provide enough information to describe the five S—P mixing amplitudes discussed in Sect. [3l 
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For TV stemming from the QCD 9 term, the proper account of chiral symmetry radically changes the 
picture. As noticed in Ref. [2U] , at leading order the 9 term generates only the isoscalar pion-nucleon T- 
odd coupling go, and thus contributes at tree level only to the I — potential [29]. A coupling of g\ form 
appears two orders down in the ChPT expansion, and the one of <?2 form is even more suppressed [25j . To 
evaluate the effects of the 9 term on observables which, like the deuteron EDM, are mostly sensitive to 
the I — I components, it is necessary to consider the TV two-nucleon potential to next-to-next-to-leading 
order in ChPT. As described in Sect. 2] this implies the consideration not only of the non-derivative 
pion-nucleon TV couplings, but also of subleading TV derivative couplings, of power-suppressed TC 
interactions (with particular care for isospin-breaking operators, which contribute to the I — 1 and 1 = 2 
potentials), and of one- loop and short-range contributions to the two-nucleon potential. When all these 
elements are considered, the potential has a much richer structure than Eq. ([§2"]) : (i) a hierarchy emerges 
between the five spin-isospin structures already present in Eq. (|92p . and (ii) momentum-dependent 
potentials appear, with the same importance as most of the momentum-independent ones. 

We first analyze the implications of our 9 results to V^ min (r). Using the chiral index v : as defined in 
Eq. (p~5|) . to keep track of the size of different pieces, the 9 contributions to Eq. (f9"2"]) are 
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where U(r), W(r), X(r), and Y(r) are defined in Eqs. ([75]), ([75]), ([77]), and ([75]), respectively. In Eq. §F5§ 
the use of the definition ([75]) to express Vq in terms of the physical pion masses introduces subleading 
corrections in the v = 1 term, which strictly speaking would use a function U(r) that only depends on a 
common pion mass, say the neutral one, Uo(r) = exp(— m^ar) / Aur . In Eqs. ([94]) (J98j we can neglect the 
pion mass difference in U(r), and use Uq{t), the error thus introduced being at higher orders in the ChPT 
power counting. Similarly, in Eq. (|93p the use of <?a and go with their m 2 corrections included accounts 
for some v = 3 corrections, while whether or not such m 2 corrections are included in Eqs. (|94p - (|98[) 
is beyond the order we consider. In Eqs. (1941) and (|98[) we replaced Sm^, the quark- mass-difference 
contribution to the nucleon mass splitting, with the physical value of the nucleon mass splitting itself, 
m n — trip , the difference again being a higher-order contribution in ChPT. 

As one can see, at order v — 3 in ChPT all the possible spin-isospin structures considered in Refs. 
35, 32 33, 34 appear. The dominant component is the isoscalar Vo 29 . In Figs. El and we plot, 
respectively, the momentum-space and configuration-space expressions for Wo- The dashed line repre- 
sents the leading-order W^, Eqs. (l64l) and ([79]) with the use of the definition ([75]) for U(r) to express 
Vq 1 ^ in terms of the physical pion masses. The dashed-double-dotted line illustrates the effect of the 
difference between the leading OPE potential computed with the function U(r) and with Uo(r). Other 
isospin-breaking corrections, which come from the nucleon mass splitting in Vq , Eqs. (|69[) and (|83[) , arc 
very small, as indicated by the long-dashed-dotted line barely distinguishable from the x-axis. At next- 
to-next-to-leading order, Vo also exhibits a medium-range component originating in TPE diagrams and 
a short-range component. The dashed-dotted line depicts the non-analytic piece of the TPE diagrams, 
Eqs. (pOj) and ([55]) . We estimate the short-range potential by assuming the coefficient C2 in Eq. ([55]) 
to be dominated by the ln/z 2 /m 2 term with fx = itin- The rationale is that there is no obvious reason 
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Figure 5: Components of the 9-term two-nucleon potential |<f|Vo, in units of gA9o/ F^m^, as a function 
of the transferred momentum \q\, in units of m^. The (blue) dashed line denotes the leading-order OPE 
contribution with physical pion masses; the (orange) dashed-double-dotted line shows the effect of the 
pion mass difference on the leading OPE contribution; the (dark green) long-dashed-dotted line accounts 
for the even smaller effect of the nucleon mass splitting; the (purple) dashed-dotted line is the non-analytic 
TPE contribution; and the (green) dotted line presents an estimate of the short-range component of the 
potential. The (red) solid line is the sum of all contributions up to next-to-next-to-leading order, except 
for the short-range component. 



to expect that such a contribution, non-analytic in m^, should get exactly canceled by m^-independent 
short-distance contributions. However, the sign cannot be guessed reliably and our choice is purely ar- 
bitrary, for illustration only. Equation (JSHJ) gives rise to the straight dotted line in Fig. [S] but Eq. ([55)1 
does not appear in Fig. [5] since it is concentrated at r — 0. The solid line in both figures is the sum of 
the long and medium range contributions to Wo- 

From Fig. [5J we can appreciate that, as expected from the ChPT power counting, the medium 
and short-range corrections to the TV potential have comparable size in the momentum range we are 
considering, and for momenta q > they noticeably affect the leading order. At momenta of order 300- 
400 MeV the medium and short-range contributions have roughly the same size as the leading potential. 
In this region, degrees of freedom which we have not explicitly included in the EFT, like the A isobar, 
become relevant, and the convergence of the perturbative expansion can be improved by extending the 
EFT to incorporate them. Isospin-breaking, long-range corrections, although of formally the same order 
as TPE and contact terms, are much smaller, at least in part because of factors of e, except at very small 
momenta where their longer range compensates. In Fig. [5] we focus our attention on the long-distance 
region, r > l/m n . At distances of up to r < 2/m,r TPE is still the dominant correction to the potential, 
but it is overcome at longer distances, r > by the long-range effects of pion mass splitting in OPE. 

It is instructive to compare our results for Vo to the corresponding TV potential obtained in a one- 
boson-exchange model. In such a model, T violation in the coupling of a rho meson to the nucleon 
generates corrections to Vq of the form [33J GH] 



y(p) (r) = _9_a^9pNN_ 9^F JL e__^ 
F£ 9-kNN 5o 47rr 



2G 




Figure 6: Components of the (9-term two- nucleoli potential Wo in units of gAgorn\/AnF% : as functions 
of the distance between the two nucleons r = |rj, in units of Curves as in Fig. [5j except that the 

short-range component of the potential is not shown. 



where g P NN is the TC rho-nucleon vector coupling and go p is an isoscalar, TV, one-derivative rho-nucleon 

coupling, defined, for example, in Ref. [34]. In the limit where the rho mass is large, m p oo, Vg^(r) 
approximates a delta function and the effect of TV rho-exchange amounts to a contribution to C 2 of the 
form 

d ( P ) = _ 2 9A9Q9 P NN9 0p F n l 
F4 g^NN go 

Since m p ~ Mqcd and there is no reason for go p F v /g to be particularly big or small, the size of rho- 
meson contribution is comparable to the power-counting expectation, C% = 0(6m%/ F%Mq CD ), with 
some suppression coming from the numerical smallness of the TC rho-nucleon vector coupling compared 
to the pion-nucleon coupling. Assuming the TV pion-nucleon and rho-nucleon couplings to have the same 
strength, gopF^Jgo = 1, and using for the rho-nucleon vector coupling the value determined in modern 
high-precision two-nucleon potentials, g P NN = 3.2 [83) . in Fig. [7] we compare the rho- meson contribution 
to Wo to the pion-mass-splitting and TPE medium-range corrections discussed above. For r > l/m w , 
the contribution of the rho meson is numerically small compared to both pion-mass-splitting and TPE 
corrections. At shorter ranges, r < l/m w , rho exchange overcomes the effect of pion mass splitting, but it 
always remains smaller than TPE. Of course we can make one-rho exchange more important by jacking 
U P SopFw/go, but we cannot compensate for the different ranges of the two contributions, m p versus 2m T . 
We see little justification for neglecting TPE in the 9 potential. 

We now turn to the other spin-isospin structures in Eq. (|92[) . which in EFT are all suppressed by 
one power of Q 2 /Mq CD with respect to the leading OPE TV potential. The function only receives 
contributions from short-range physics. Again, in a one-boson-exchange scenario, contributions of exactly 
the size of C\ come from eta and omega exchanges [3H [33J EH] , 

£i(J7,w) _ 2^450 f 9r]NN §0riFn 1 ffcjjVAf gOcjfir 1 \ (101) 

1 F% \girNN go g^NN go ml J ' 

where m v (m u ) is the eta (omega) mass, g^NN {9uiNn) is the TC eta-nucleon axial (omega-nucleon vector) 
coupling, and g^ (gow) is an isoscalar, TV no-derivative eta-nucleon (one-derivative omega-nucleon) 
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Figure 7: Comparison between one-rho-exchange and EFT contributions to the magnitude of the 8 
potential |VVo| in units of gAgarn^/AirF^, as functions of the distance r, in units of l/m T . The rho- 
exchange contribution is depicted as a (blue) long-dashcd-dottcd line, while TPE and pion mass splitting 
in OPE are as in Fig. El 



coupling. The eta- and omega-meson contributions are comparable to the power-counting expectation 
C\ = Q(6m^l ' F%Mq CD ). For the eta meson, the enhancement due to the relatively light mass is offset 
by the smallness of eta-nucleon TC coupling, g V NN = 2.24 [84]. The ratio g^NN / g-KNN is instead close to 
one (83], and, therefore, we have no reason to expect the omega- meson contribution to C\ to differ much 
from the power-counting estimate. 

(3) (3) (2+3) 

In contrast, V{ ,U\ , and V 2 sprout entirely from OPE. The TV, isospin-breaking coupling g\ 
contributes equally to Vi and Ui, as expected [30] from the identification at the Lagrangian level, cf. 
Eqs. (|24|) and (|90|) . However, we expect a comparable long-range piece in U\ — Vi, which stems from 
the combination of the isospin-violating vertex fii and the TV vertex go. As discussed in Ref. |25) . 
strong-dynamics contributions to the coefficients of these 1=1 potentials are in principle determined 
by measurement of TC, isospin-breaking observables. For example gi/go could be extracted from a 
detailed analysis of isospin-breaking effects in pion-nucleon scattering. At present, however, even the 
very sophisticated, state-of-the art analysis of Ref. [85] stops one order shy of the accuracy required 
for such extraction. Similarly, the ratio fii/gA affects isospin violation in nucleon-nucleon scattering, 
but at present phase-shift analyses of two-nucleon data can only provide a bound on Pi, which is in 
accordance with the power-counting expectation [51] [52] . In the absence of better constraints on the 
parameters in Eqs. ([96]) and ([97]) . the ratios Ao an< ^ ^i^Ao^ can om y be estimated by power 
counting, as 0(em^/MQ CD ) ~ 1%. As for the last component of the phenomenological potential, V^ 3 ' 
originates entirely from the isospin-violating corrections to the pion and nucleon masses, and it is also 
relatively small. Note that to this order this component has nothing to do with the coupling g~2 of a 
phenomenological Lagrangian: because of the isoscalar character of the 8 term, Eq. ([2~T]) . gi arises in 
EFT only at higher order. 

In one-boson-exchange models the I — 1,2 potentials are assumed to arise from pion, eta, rho, and 
omega isovector and tensor TV couplings to the nucleon [32] [33l [34]. In the ChPT power counting, 
short-range contributions to these potentials are suppressed with respect to the long-range pieces, again 
because of the isoscalar character of the 8 term. (Of course, because of the factors e in the long-range 
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contributions of this order, short-range terms might not be entirely negligible.) This is consistent with 
the argument that the dominant meson-exchange contributions are from the pion and the eta [32] . 

There are, therefore, a few points of contact between the local part of our v < 3 potential and the 
phenomenological potential Vp min (r) (|92|) . However, as we have seen in Sects. H]and[5j at this order EFT 
yields also momentum-dependent interactions, which in coordinate space appear as non-local potentials 
and corrections that account for CM motion of the nucleon pair. They can be found in the relativistic 
and isospin-breaking corrections to OPE in Eqs. (|79[) . (I5TT) . and (|82l) . 

At v = 3, the 9 two-nucleon potential contains in the CM frame, P = 0, four spin-isospin structures 
that are momentum-dependent, 



-\ (V^(r)) x a< 2 )) ■ p r + (c?W x o^)™ Uj™V r U{r) - ^"V 2 .^)) p\ 



+B ^^^^ + ^y^ u{r)h (102) 

where p r = — iV r denotes the quantum-mechanical relative momentum operator. 

The structure of the momentum-dependent TV potentials was considered previously in Ref. [35] , where 
all possible Hermitian operators were constructed, which violate time-reversal and parity, and contain up 
to one power of momentum p r . The momentum-dependent TV potential was parameterized with eleven 
unknown functions di(r), i — 1,2..., 11. The first term in Eq. (|102[) is quadratic in the momentum 
operator and was not considered in Ref. [35] . The second and third spin-isospin structures correspond, 
respectively, to the isoscalar functions d,2(r) and de(r). For TV from the QCD 9 term, these two functions 
are therefore dominated by pion-exchange, and their coefficients are fixed by Lorentz invariance and do 
not contain any new TV parameter. Isospin-breaking effects in the strong interaction give rise to the last 
term in Eq. (I102[) . which is proportional to the nucleon mass difference, and it is the first contribution 
of the 9 term to dio(r). Once again, dio is dominated by OPE diagrams, and the only TV parameter 
intervening is go. We find that, at order v — 3 in ChPT, the other isospin-conserving (the isoscalar d\ 
and d§) and isospin-breaking (d^, c?4, c?7, rfg, dg and du) functions do not receive contributions from 9. 

In order to get a sense of the importance of the momentum-dependent contributions, we consider the 
effect of the relativistic correction in Eq. (I102[) that is quadratic in p r . In Fig. |S]we compare it (long- 
dashed-dotted line) to leading OPE (dashed line), medium-range TPE (dashed-dotted line), and pion 
mass-splitting corrections (double-dotted-dashed line), all applied to a simple bound-state wave function 
with the scale present in the 1 S'o channel, a s — —23.714 fm: 

= exp(-r/ fl .) 
r 

In this qualitative example the relativistic correction cannot be neglected with respect to the other v = 3 
corrections, and we take Fig. [8] as an indication that also in actual calculations of TV observables it 
would not be advisable to drop the potential in Eq. (|102l) . when the 9 potential is needed to next-to- 
next-to-leading-order accuracy. 

Finally, to the same order we find contributions proportional to the CM momentum of the nucleon 
pair, 



+ (a^-a^y 2(V r *7(r))p 2 + P (v r C/(r)) • P + i(P • V r ) 2 (v>t/(r)) j 



(104) 
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Figure 8: Comparison between a relativistic correction to OPE and local components of the 9 two- 
nucleon potential Wo applied to an illustrative bound-state wave function ip, in units of gAgorn^ / AnF% , 
as functions of the distance r, in units of l/m v . The (dark green) long-dashed-dotted line represents the 
term in the potential that is quadratic in momentum. Other curves are as in Fig. [5] 

where P = —iVx- Although the operators in Eq. fll04[) vanish in the two-nucleon CM frame and are 
not important for the study of T violation in nucleon-nucleon scattering, they impact observables like 
TV electromagnetic form factors of the deuteron, where the recoil against the photon changes the CM 
momentum of the nucleon pair, and they have to be considered in nuclear systems with A > 2. An 
example of the effects of recoil on TC deuteron processes can be found in Compton scattering 86 . 

At leading order, the dimension-6 sources of TV only contribute to VJr. m ; n (r), the two-nucleon potential 
with minimal number of derivatives, in Eq. (|92[) . Combining all these sources, 
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where 

• for qCEDM, f m in = —1 and only the <?o,i terms apply; 

• for TVCI sources, v m \ n = — 1 and only the go,i an d Ci,2 terms apply; and 

• for qEDM, f m i n = 2 and only the go, 1,2 and c?o,i terms apply. 
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In the case of the qCEDM, the potential is thus dominated by OPE, with shorter- or longer-range 
contributions expected to be small. For the gCEDM and the TV FQ operators, short-range effects are 
also leading; as for 9, these can be parametrized by rho (see Eq. (llOOp ). and eta and/or omega exchange 
(see Eq. f|101[) ^ . Only for qEDM all components of Vr tin ia(^) appear in leading order thanks to both the 
most general non- derivative pion-nucleon coupling structure, and to the long-distance potential stemming 
from the nucleon EDM. The latter cannot be well approximated by heavy-meson exchange. 

7 Conclusions 

The power-counting scheme of EFTs allows us to organize the contributions to the potential in powers 
of Mq^ d . The TV PV potential has an ordering that has similarities with the TC P V potential [56] : for 
all TV sources of dimension up to 6, the leading potential contains one-pion exchange. For T violation 
from the 9 term and from the quark chromo-EDM, the relative importance of two-pion exchange and 
shorter-range interactions follows the TC PV case closely. However, the situation is different for the quark 
EDM, the gluon chromo-EDM and TV four-quark operators, where contact interactions or one-photon 
exchange are relatively more important. 

For T violation from the 9 term, at leading order, 0(Q/Mqcd), we find only the well-known OPE 
from the 1 = pion-nucleon TV coupling [29]. The OPE from the 1=1 pion-nucleon TV coupling is 
suppressed by two orders in the expansion parameter and is of 0(Q 3 /Mq CD ). Since the 1 = OPE is 
suppressed in nuclei, higher orders in the potential could be important. We have thus also examined the 
corrections in the next two orders, which are up to 0(Q 2 /Mq CD ) relative to leading. We have found 
that the potential is purely two-body, and: 

• At the longest, one-pion range, there are more general vertex corrections than usually assumed. We 
employed the results of Ref. |25] where the TV pion-nucleon vertex was examined to this order. 
In addition to the qualitatively different 1=1 pion-nucleon TV coupling, there are corrections to 
the local potential stemming from isospin breaking in the pion and nucleon masses, and in the TC 
pion-nucleon coupling. There are also recoil (oc 1/m^) and relativistic (oc l/m 2 N ) corrections to 
leading OPE, which make the potential non-local and dependent on the total momentum of the 
nucleon pair. 

• At this order, we find, additionally, two-pion exchange from the 1 = TV coupling. The non- 
analytic, medium-range part of the TPE potential is independent of the choice of fields and regula- 
tors. Like the leading OPE potential, this part of the TPE potential has as only (so-far) unknown 
quantity the 7 = TV pion-nucleon coupling. Under reasonable assumptions about the strengths 
of TV couplings, this potential is stronger, and has a different radial dependence, than phcnomeno- 
logical one-meson-exchange potentials. The main effect of TPE is to modify the potential in the 
same channels as the leading OPE potential. 

• The short-range part of the TPE potential, on the other hand, cannot be separated from contact 
interactions, the most general form of which we also write at this order. They are two of the terms 
given in the literature [34) . In the context of a theory without pions, this implies a contribution to 
only two of five possible S-P transitions. When pions are included explicitly, the short-range terms 
are expected to be of the same size as TPE, and thus their strengths depend on the renormalization 
scale. They subsume short-range dynamics that includes the effects of heavier mesons, but whether 
such effects are sufficient to saturate them is unknown. 

The structure of the resulting potential is therefore significantly different than the phenomenological 
potential used in the literature, due to the specific way in which the 9 term breaks chiral symmetry. 
If consideration of short-range dynamics or J = 1 OPE is necessary, one should also include the TPE 
potential and OPE corrections calculated here. 
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For the dimension-6 sources, the 1=1 pion-nucleon coupling appears already in leading order, 
0{M QCD /Q) for qCEDM and chiral-invariant sources (gCEDM and TV FQ) and 0(Q 2 /M^ CD ) for 
qEDM. The structure of the TV potential is thus different from that of 9, and depends on the source: 

• For qCEDM, it contains only OPE with the well-known 7 = 0,1 non-derivative pion-nucleon cou- 
plings. 

• For CI sources, there are, additionally, two contact interactions of the same type as found at next- 
to-next-to-leading order for 9. 

• For qEDM, OPE from the three {I = 0, 1, 2) non-derivative pion-nucleon couplings is accompanied 
by one-photon exchange from the nucleon EDM. 

Therefore, for all sources we have considered (9, qEDM, qCEDM, gCEDM, and TV FQ), the TV 
nuclear potential presents more structure (spin, isospin, and/or distance profile) than assumed in the usual 
phenomenological approach. EFT offers a framework where the calculation of nuclear TV observables 
can be carried out in a model independent way and characterize the low-energy manifestations of possible 
TV sources. 
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Appendix: Fourier Transformation to Configuration Space 

In general a potential obtained in EFT depends not only on the transferred momentum q but also on K, 
and the CM momentum P, V(q, K, P). The Fourier transform of such a potential is defined as 

V(f,f'Xx') =J^f**LJ ^e^<*-*^We-i*Wvtoit,?), (HO) 

where, if X\ and xi are the positions of the incoming nucleons and x[ and x' 2 the positions of the outgoing 
nucleons, the relative coordinates are r = X1 — X2 and f ' = x{ — x*2, while the CM position of the incoming 
and outgoing pairs are 2X = X\ + X2 and 2X 1 = x[ + x 2 - The potential in Eq. (jllOl) has to be used in a 
two-nucleon Schrodinger equation of the form 

i^(f',X') = - ^L + ^j iP(r',X') + J d 3 fj d 3 XV(r,r',X,X')iP(r,X). (Ill) 

For potentials that, like the ones in Sect. 21 are polynomials in K and P, 

V{q,K,P)^K m P n f(q), (112) 

V(r, r ', X, X 1 ) assumes the form 

V(r,?',X,X')oc (V^(X-X')) (v?5V\f-n)f(p^-), (113) 

where 

f ^= I 7?T3 e ^/(9)- (114) 

J ( Z7T) 
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Plugging Eq. (|113p in Eq. and integrating by parts, the derivatives acting on the delta functions 

can be turned into derivatives acting on / and on the wave function ip(r, X). The integrals in Eq. (jllip 
then become trivial, and the Schrodinger equation assumes the form 

i^(r',X') = - (p^- + ^ j iP(r',X') + V(r>,V r ,,V x ,)iP(r>,X>), (115) 
at \ Am A? rriM I 

where the two potentials in Eqs. (|llip and (| 1 15[) are related by integrations by parts. For a potential of 
the form (|112p , schematically we would have 

V(r\V r ,,V x ,) ex (-) n (-r [^Y L , -- V£„ (116) 

where the indices ij, . . . , i m are appropriately contracted. 

In order to obtain the potential in configuration space for functions that diverge as the momentum 
transfer |g| goes to infinity, one has to define a regularization scheme. Here, following Ref. [82], we find it 
convenient to extend the definition of the Fourier transform (| 1 14[) to a space-time of d = n+l dimensions: 

V n {r) = J -^e-^V(q). (117) 

The amplitude V(q) is the expression in momentum space of corresponding loop contributions right after 
the d-dimensional integration over loop momenta is performed, but before setting d — 4 or performing 
the integration over Feynman parameters. Writing 

d n q = q^dqil - cos 2 6)^ dcos0dfi„_i, (118) 
the angular integrations are evaluated with the aid of the formulas 
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. dcos9 (1- cos 2 0)^ e -i«rco S e = (27r )f (qrf-^J^qr), (120) 

1 {—) J-l 

where q = \q\, r = |r |, and J n (x) denotes a Bessel function of the first kind. For momentum integrals, a 
useful relation is [87] 

n Ja.-x(qr) /r\ 1JL r 1 - 1 8™-% 
dqq* 2 n V-n = UJ FT^Y^-nW, (121) 

(g 2 + /? 2 )— V27 r(^) 

where /? is a constant and if „ (x) is the modified Bessel function of the second kind. 
For example, in the case of the triangle diagrams discussed in Sect. |U 

v Ai q) - ■ - (2) - - (2) ) • rr (^) jf 1 & K + ?Mi - x)] ^ • (122) 

For m = 3 and /3 2 = to 2 /x(1 — x), the result in Eq. (|121l) allows one to cancel the divergent factor of 
F((3 — n)/2) in Eq. (|122l) and get an expression that is finite for r ^ 0. Now we can set d = 4 and 
with the aid of the properties of modified Bessel functions |87j we can write the potential in configuration 
space as 



1 f 1 



27rr 3 



dx(l + Br) e~ 0r 



(123) 
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The contributions from box and crossed diagrams can be obtained in a similar fashion, leading to the 
result in Eq. ([55]). 

Alternatively, we can isolate the short-range, divergent part of the interaction with integration by 
parts. In Eq. ()122|) . for instance, we then obtain 



.gA90 (4tt^ 2 ) 2" , (2)/-(l) -.(2)\ - 



3 — n 



5 — n 



dx- 



q 2 x{\ - 2x) 



\m% + q 2 x(l - x)]' 



(124) 



The first, divergent piece in Eq. (| 124[) is a contribution to a delta- function potential. Applying the 
rf-dimensional Fourier transform to Eq. (|124[) and taking the d —> 4 limit we find 



V A {r) 



9a9o 



F 2 (2nF 7r ) 2 



(1). T (2) (#(l)_ 3 l2)y$ 



1 - 2x 



47rr J 1 — x 



-/3r 



(125) 



where L is given in Eq. (|57p . Proceeding in this way also for box and crossed terms, we find that Fourier 
transform of the TPE potential can be expressed as 



1 



F* (2^)2 

rl 

dx 
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-(1) . T (2)(£(l) _^(2)).V J-,5( 3 )(f) 



(3. 9 ; 



l)(L + lnAr 



2 5 



1 

47rr 



2x(l - x) 



1 - 2a; 



1-x 



2 e -Pr 



(126) 



The piece proportional to the delta function can then be absorbed in a redefinition of Ci very similar to 
Eq. (|58p. the only difference residing in the finite pieces. Integrating by parts, it can be explicitly verified 
that the non-analytic piece of the expression (|126j) gives the medium-range potential in the form of Eq. 

(113. 

As a further check of our results, we computed the Fourier transform of the triangle diagrams with 
a Gaussian regulator exp(— q 2 /A 2 ), for different values of the cutoff A. The calculation was performed 
numerically with Mathematica [88j and we focused on the region r > l/m^. For A ~ m p , the result we 
get is still quite different from the Fourier transform obtained in dimensional regularization, but as we 
increase the cutoff to 1-2 GeV, it approximates Eq. (|123[) better and better. 

As pointed out in Ref. [82] , in the d-dimensional Fourier-transform procedure the infinities are "regu- 
larized" away because the nucleon distance is kept finite. The ultraviolet divergences and the regulator 
dependence are now hidden in the singular behavior (~ 1/r 4 ) of the potential for small r, which forces 
the reintroduction of a regulator in the calculation of matrix elements of V(r). If the chosen regulator 
were dimensional regularization, then the d 4 limit of Eq. (|121j) . which leads to the 1/r 3 singularity 
in Eq. (I123p . must be taken in the sense of generalized functions; the singularity is then encoded by a 
delta function, proportional to the divergent factor 2/{d— 4) , with a plus distribution remaining [89 . 
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